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ABSTRACT 

We investigate the clustering properties of 2,600 Lyman Break Galaxies (LBGs) at z = 3.5 — 5.2 in 
two large blank fields, the Subaru Deep Field and the Subaru/XMM Deep Field (600arcmin 2 each). 
The angular correlation functions of these LBGs show a clear clustering at both z ~ 4 and 5. The 
correlation lengths are r = 4.1+°:2 and S^tx^ioo M P C ( r o = and ^^-x^Kao M P C ) for a11 tne 

detected LBGs (for L > L* LBGs) at z ~ 4 and 5, respectively. These correlation lengths correspond 
to galaxy-dark matter biases of b g = 2.91°^ and 4.6t°;2 ( & s = 3 - 5 io.7 and 4 - 6 iiJ)i f° r a11 tne detected 
LBGs (for L > L* LBGs) at z ~ 4 and 5, respectively. These results, combined with estimates for z ~ 3 
LBGs in the literature, show that the correlation length of L > L* LBGs is almost constant, ~ 5^ioo 
Mpc, over z ~ 3 — 5, while the bias monotonically increases with redshift at z > 3. We also find that for 
LBGs at z ~ 4 the clustering amplitude increases with UV-continuum luminosity and with the amount 
of dust extinction. We estimate the mass of dark halos hosting various kinds of high-z galaxies including 
LBGs with the analytic model given by Sheth & Tormen (1999). We find that the typical mass of dark 
halos hosting L > L* LBGs is about 1 x lO 12 /if o 1 M over z ~ 3 — 5, which is comparable to that of the 
Milky Way Galaxy. A single dark halo with ~ lO 12 /if o 1 M is found to host 0.1-0.3 LBG on average but 
host about four if-band selected galaxies. 

Subject headings: cosmology: observations — cosmology: early universe — cosmology: large-scale 
structure of universe — galaxies: high-redshift — galaxies: evolution 



1. INTRODUCTION 

Formation history of galaxies is basically understood by 
two fundamental evolutionary processes, i.e., production of 
stars and accumulation of dark matter. Electromagnetic 
radiation from galaxies gives information about the status 
of baryonic matter, i.e., stars, gas, and dust. Although 
dark matter should play a key role in galaxy formation, it 
cannot be directly detected using electromagnetic waves. 
Clustering properties of galaxies are closely related to the 
distribution and amount of underlying dark matter (e.g., 
Peacock et al. 2001; Percival et al. 2001; Verde et al. 2002; 
Percival et al. 2002; Szalay et al. 2003). Thus, measuring 
clustering properties of galaxies is very useful for under- 
standing the dark side of galaxy properties (different from 
stellar properties). 

Clustering properties in the present-day universe have 
been well investigated using recent wide-field surveys such 



as 2dFGRS (Colless et al. 2001) and SDSS (York et al. 
2000). Clustering properties of galaxies are derived from 
spatial or angular two-point correlation functions (CFs), 
and are compared with those of the dark matter predicted 
by Cold Dark Matter (CDM) models (Peacock et al. 2001; 
Percival et al. 2001; Verde et al. 2002; Percival et al. 2002; 
Scranton et al. 2002; Connolly et al. 2002; Dodelson et al. 
2002; Tegmark et al. 2002; Zehavi et al. 2002; Szalay et al. 
2003). An analysis of the power spectrum from 2dFGRS 
data shows that the distribution of galaxies agrees with 
that of dark matter (Verde et al. 2002). In other words, 
optically selected galaxies trace the underlying mass dis- 
tribution at z = 0. 

Clustering properties of galaxies are now investigated 
at higher redshifts up to z ~ 4. Daddi et al. (2000) have 
reported a strong angular correlation for EROs at z ~ 1, 
and estimated the correlation length of the spatial CF to 
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be r = 12 ± 3h^ Mpc, which is comparable to that of 
present-day ellipticals (Daddi et al. 2001). McCarthy et 
al. (2001) found similarly large, but slightly smaller, corre- 
lation lengths (r ~ 9 — lO/ifoo Mp c ) f° r their red galaxies 
at z ~ 1. Miyazaki et al. (2003) found strong clustering for 
both old galaxies and dusty star-forming galaxies beyond 
z ~ 1. Giavalisco et al. (1998) have studied CFs of Lyman 
Break Galaxies (LBGs) on the basis of a large sample for 
bright (R < 25.5) LBGs lying around z = 3 selected by 
U n GR colors. They have found that the spatial distribu- 
tion of z = 3 LBGs is strongly biased relative to the dark 
matter distribution predicted by CDM models with a lin- 
ear bias of 4.5 for an Einstcin-de Sitter cosmology (see also 
Adclberger et al. 1998). In a subsequent paper, Giavalisco 
& Dickinson (2001) have found the clustering amplitude 
of z ~ 3 LBGs to depend on their rest-frame ultraviolet 
luminosity, with fainter galaxies less strongly clustered, 
which is a similar property to that found in the present- 
day universe. Arnouts et al. (1999, 2002) have measured 
the correlation length of faint galaxies in the HDF-N and 
HDF-S over the redshift range < z < 4, though er- 
rors in their measurements are large due to small number 
statistics. Recent observational results show that LBGs 
at z ~ 4 and Lyman a emitters (LAEs) at z — 4.9 are 
strongly biased against underlying dark matter, and that 
the clustering amplitudes of LAEs are segregated with re- 
spect to Lya luminosity (Ouchi et al. 2001; Ouchi et al. 
2003a, hereafter SDS II). 

Narrow-band searches of LAEs on targeted fields have 
found strong clustering of LAEs around quasars (Cam- 
pos et al. 1999; M0ller & Fynbo 2001; Stiavelli et al. 
2001; Fynbo et al. 2003), radio galaxies (Venemans et 
al. 2002), high redshift clusters (Palunas, Francis, & 
Woodgate 2000), and overdensities in the LBG distribu- 
tion (Steidel et al. 2000). 11 Recent deep and wide-field 
LAE searches with Subaru/Suprime-Cam have revealed a 
large-scale distribution of LAEs (SDS II), and Shimasaku 
et al. (2003, hereafter SDS IV) have found a filamentary 
large-scale structure made of LAEs at z = 4.9 whose width 
and length are 20h^Q Mpc and > 50h^Q Mpc, respectively. 
The size of this high-z large-scale structure is compara- 
ble to those of present-day large-scale structures found by 
Gcller & Huchra (1989). The existence of such a high-z 
large-scale structure also indicates that galaxies at z ~ 5 
would be highly biased against underlying dark matter. 
However, LAE searches by narrow-band imaging observe 
a thin slice of high-z universes, resulting in a surveyed vol- 
ume of as small as about an order of 10 4 ft.f 3 Mpc 3 . Such 
surveys may suffer from a strong field variance. On the 
other hand, survey volumes for LBGs are generally about 
20-50 times larger than those for LAEs. Thus, measuring 
clustering properties of LBGs is highly desirable to quan- 
tify the average distribution of high-z galaxies. 

Motivated by this, we carried out deep and wide-field 
imaging for two blank fields, the Subaru Deep Field 
(SDF: 13 h 24 m 21.4 s ,+27°29'23"[J2000]; Maihara et al. 
2001, hereafter SDS I) and the Subaru/XMM Deep sur- 
vey Field (SXDF: 2 /l 18 m 00 s ,-5°00'00"[J2000]; Sekiguchi 
et al. 2003 in preparation; see also Ouchi et al. 2001), and 



made samples of 2,600 LBGs at z = 3.5 — 5.2 distributed in 
these two fields (600 arcmin 2 each). Details of the samples 
are described in Ouchi et al. (2003b, hereafter SDS V). 

In this paper, we investigate clustering properties of 
LBGs at z = 4 and 5 using these large LBG samples. 
We then compare the observational results with predic- 
tions of CDM models as follows. Applying analytic mod- 
els for the spatial clustering of dark matter to the obser- 
vational results, we estimate the mass of dark halos which 
host galaxies from the bias of galaxy-dark matter distribu- 
tion (Bullock, Wechsler, & Somerville 2002; Moustakas & 
Somervillc 2002). Since the dark- halo mass is a fundamen- 
tal quantity of galaxies, we next compare various high-z 
galaxies detected in various wavelengths in terms of halo 
mass. In reality, a large fraction of high-z galaxies de- 
tected to date are either LBGs, LAEs, EROs, or SCUBA 
sources. These high-z galaxy samples are biased toward 
UV-bright star- forming galaxies (LBGs), strong- Lya emis- 
sion galaxies (LAEs), old passive galaxies (EROs), and 
dusty starburst galaxies (SCUBA sources). It is essential 
to understand the nature of the high-z galaxy population 
as a whole by combining these biased samples in order to 
propose a scenario for the evolution of galaxies at high 
rcdshifts. Finally, we examine the relation between high-z 
galaxies and present-day galaxies on the basis of dark-halo 
mass of descendants predicted by the CDM model (Mous- 
takas & Somerville 2002). 

The outline of this paper is as follows. In section 2, we 
describe photometric samples of LBGs at z = 3.5 — 5.2 
from the data in the SDF and the SXDF, which are ob- 
tained by SDS V. These are the largest samples of LBGs 
at z > 4 in contiguous areas obtained to date, and thus 
enable detailed studies of clustering properties of galaxies 
at the highest redshifts so far. Section 3 presents obser- 
vational results. We investigate the spatial distributions 
of LBGs to obtain the observational results. We derive 
angular correlation functions (CFs) of LBGs, and calcu- 
late correlation lengths of spatial clustering. We estimate 
the galaxy-dark matter bias of LBGs at z ~ 4 and 5, and 
compare them with those of low-z galaxies and LBGs at 
z = 3. We examine the dependence of clustering ampli- 
tude on luminosity and color for LBGs and LAEs. Section 
4 presents theoretical implications. We compare our ob- 
servational results on the CF and the luminosity function 
obtained in SDS V with an analytic CDM model to inves- 
tigate properties of dark halos hosting high-z galaxies and 
the relation between galaxies (luminous matter) and dark 
halos. In Section 4, we also infer the mass of descendants 
of these high-z galaxies. In section 5, we propose a uni- 
fied view for a variety of high-z galaxies and discuss the 
formation history of galaxies from z = 5 to z = 0, based 
on the results obtained in this paper and SDS V. Section 
6 presents a summary of this paper. 

Throughout this paper, magnitudes are in the AB sys- 
tem (Oke 1974; Fukugita, Shimasaku, & Ichikawa 1995). 
The values for the cosmological parameters adopted in this 
paper are: (0 TO , n, eg) = (0.3,0.7,1.0,0.9). These val- 
ues are the same as those obtained from the latest CMB 
observations (Spergel et al. 2003). 12 We adopt n = 1 



11 Searches of LBGs in targeted fields are also made by, e.g., Bouche & Lowenthal (2003) who find a significant clustering of LBGs around a 
damped Ly a absorber at z = 3. 

12 Spergel et al. (2003) show that the best-fit parameters arc n = 0.93 ± 0.03 and ag = 0.84 ± 0.04. However, the differences between these 
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and as = 0.9, since these values were used in many of the 
theoretical and observational studies to date. We express 
physical quantities using h 70 , where /170 is the Hubble con- 
stant in units of 70 km s _1 Mpc -1 . The exception is the 
correlation length, r . We express r using h wo , where 
is the Hubble constant in units of 100 km s _1 Mpc -1 , 
since most of the previous studies use ft-ioo to express r . 

2. OBSERVATIONS AND GALAXY SAMPLES 

We use three Lyman Break Galaxy (LBG) samples pre- 
sented in Ouchi ct al. (2003b), i.e., SDS V. Details of the 
samples are shown in SDS V, and we give a brief summary 
of the data and the samples in the following. In sections 
3-4, we describe the clustering properties of a sample of 
Lyman a emitters (LAEs) at z = 4.86 ± 0.03 found in the 
SDF to compare them with those of LBGs. Details of the 
LAE sample and their clustering properties (e.g., distribu- 
tion and angular correlation function etc.) arc shown in 
SDS II. 

2.1. Observations 

During the commissioning runs of Suprime-Cam from 
November 2000 to November 2001, we carried out 
multi-band, deep and wide-field optical imaging in two 
blank fields. One is the Subaru Deep Field (SDF: 
13 /l 24 m 21.4 s ,+27°29'23"[J2000]) covering a 616 arcmin 2 
area whose central region has very deep J and K' im- 
ages (SDS I). The other is the Subaru/XMM Deep 
Field (SXDF: 2 /l 18 m 00 s ,-5°12'00"[J2000]) covering a 653 
arcmin 2 area. We took multi-color images with the B-, V- 
, R-, i'-, and z'-band filters for both fields. The exposure 
time ranges from 81 to 210 (40 to 177) minutes among the 
bandpasses in the SDF (SXDF). The limiting magnitude 
and the seeing size of the reduced i' image is i' = 26.9 
and 0".9 (i' = 26.2 and 1".0) for the SDF (SXDF), re- 
spectively. We limit the object catalogs to i' < 26.5 and 
z' < 26.0 {i 1 < 26.0 and z' < 25.5) for the SDF (SXDF), 
in order to provide a reasonable level of photometric com- 
pleteness. Our z'-limited and z'-limited catalogs contain 
respectively 45,923 (39,301) and 37,486 (34,024) objects 
for the SDF (SXDF). 

2.2. Definition of BRi, Viz, and Riz-Lyman Break 
Galaxies 

We make three photometric samples of LBGs by the fol- 
lowing two color selections. The first sample is for BRi- 
LBGs at z ~ 4. Their Lyman break enters into the B band 
and thus they are identified by red B — R and blue R — i' 
colors. Similarly, the second is for Viz-LBGs at z ~ 5 se- 
lected by V— i' and i' — z', and the third is for i?zz-LBGs at 
z ~ 5 selected by R — i' and i' — z'. Wc define the selection 
criteria for these LBGs so that the completeness be suffi- 
ciently high and the contamination be negligibly small. In 
order to determine the selection criteria in the two-color 
diagrams, we use the best-fit SEDs of galaxies given in 
Furusawa et al.'s (2000) HDF-N photometric redshift cat- 
alog. Wc generate artificial galaxies that mimic the colors 
of the HDF-N galaxies, and distribute them randomly on 
our original images after adding Poisson noises. Then, we 
detect these simulated objects and measure their bright- 
ness in the same manner as for our original images. We 
iterate this process 100 times, and generate probability 

values and those we adopt here are less than 10%. 



maps of low-z interlopers and LBGs in the two-color dia- 
grams. Based on the probability maps, we determine the 
selection criteria for LBGs which give small contaminants 
and keep a sufficiently high completeness. The selection 
criteria adopted are 

B -R> 1.2, R-i' < 0.7, B - R> 1.6(7? - i') + 1.9 (BRi - LBGs) (1) 
V -i' > 1.2, i' - z' < 0.7, V -i' > 1.8(1' - z') + 1.7 (Viz - LBGs) (2) 
R-i' > 1.2, i' - z' < 0.7, R-i' > 1.0(i' - z') + 1.0 (Riz - LBGs) (3) 

We apply these selection criteria to our photometric 
catalogs, and find 1,438 (732), 246(34), and 68 (38) ob- 
jects for BRi-LBGs, Viz-LBGs, and ffiz-LBGs in the 
SDF (SXDF). The y-band image of the SXDF is not deep 
enough to produce a Viz-LBG sample with a small num- 
ber of contaminants by eq.(2). In order to avoid a high 
contamination rate, we adopt 

V — i' > 1.2, i'-z' < 0.7, V-i' > 1.8(i' - z') + 2.3 (4) 

as the Viz-LBG selection criteria for the SXDF. Tabic 1 
summarizes the LBGs found in our data. 

We estimate the redshift distribution, completeness, and 
contamination of the LBG samples using the probability 
maps obtained from the simulations. The redshift ranges 
are found to be z = 4.0 ± 0.5 for BRi-LBGs, z = 4.7 ± 0.5 
for Viz-LBGs, and z = 4.9 ± 0.3 for Riz-LBG. Figure 1 
shows the number-weighted redshift distributions. The ra- 
tio of contaminants to sample objects in number is 1% for 
BRi-LBGs, 26% for Viz-LBGs, and 40% for Riz-LBGs in 
the SDF. The completeness of each sample is 40-60% in the 
most efficient redshift range. The small ratio of contami- 
nants in BRi-LBGs is due to the fact that the BRi-LBG 
sample includes many faint LBGs beyond M* , which are 
much more numerous than foreground contaminants. 

2.3. Spectroscopic Follow-up Observations 

In order to examine how well these selection criteria 
isolate true LBGs, we use 85 spectroscopically identified 
objects at < z < 5 in our fields obtained by Kashikawa 
et al. (2003, hereafter SDS III), SDS IV, and SDS V. We 
find that 6 out of the 85 objects are LBGs at z > 3.5, 
and that no spectroscopically identified low-z interloper is 
included in our LBG samples. Thus, our LBG samples are 
thought to be reliable. 

3. SPATIAL DISTRIBUTION 

3.1. Sky Distribution 

Figures 2-5 show the sky distributions of BRi-LBGs, 
Viz-LBGs, and Riz-LBGs detected in the SDF and the 
SXDF. We plot all (z) ~ 5 objects, i.e., Viz-LBGs and 
i?iz-LBGs, in a single frame; Figure 3 is for the SDF and 
Figure 5 is for the SXDF. In addition to these LBGs, we 
also plot the positions of 87 Lyman a emitters (LAEs) at 
z = 4.86 ± 0.03 found by narrow-band observations (SDS 
II) in Figure 3. As found in Figures 3 and 5, some of the 
objects are selected by both Viz- and Riz-LBG (and/or 
LAE at z = 4.86) selection criteria. This is because these 
criteria select objects in similar redshift ranges: Viz-LBGs 
at z = 4.7 ± 0.5, Riz-LBGs at z = 4.9 ± 0.3, and LAEs 
at z = 4.86 ± 0.03. Table 2 shows the numbers of ob- 
jects identified by two selection criteria. These numbers 
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apparently look quite small. We examine whether these 
numbers are reasonable. The number of objects selected 
by two criteria is predicted by: 



A EST * (AT[SEL1] - AT C [SEL1]) x x <7[SEL2], 

(5) 

where iV[SELl] is the total number of objects selected by 
the first criterion, A C [SEL1] is the total number of con- 
taminants in the objects selected by the first criterion, 
Az[SEL2] and Az[SELl] are the approximate redshift 
ranges of the second and the first criteria, and C[SEL2] 
is the number-weighted mean completeness of the second 
criterion. For this calculation, we use the contamination, 
completeness, and redshift distribution derived in SDS V 
for LBGs and those derived in SDS II for LAEs. See SDSV 
for details. The estimated numbers of objects identified 
by two selection criteria are given in the fourth column 
of Table 2, and are consistent with the actual observed 
numbers within Poisson errors. This coincidence not only 
explains the relatively small numbers of objects identified 
by two selection criteria, but also ensures the accuracy of 
the contamination, completeness, and redshift distribution 
obtained by SDS V and SDS II. 

3.2. Tests for Detection Inhomogeneity 

Figures 2-5 show somewhat inhomogencous distribu- 
tions of LBGs. The surface density of LAEs has a large 
gradient over the whole image. Prior to investigating their 
clustering properties, we examine whether or not these in- 
homogeneities come from spatial inhomogeneities of the 
detection efficiency and/or photometric accuracy in the 
images. We find that the photometric zero points are ac- 
curate within 0.1 mag over the whole field of view for any 
bandpass, since PSF-like objects are found to make a sin- 
gle sharp stellar locus in any two-color plane as shown in 
SDS V. 

Then we examine whether there are spatial differences in 
the source detection efficiency by (1) calculating the num- 
ber densities of all detected objects in small (10' x 10') 
areas covering the survey region, (2) measuring the limit- 
ing magnitudes in 2700 small (40" x 40") areas for each of 
the B, V, R, i', z', and NB711 images, and (3) estimating 
the detection completeness of LBGs and LAEs from Monte 
Carlo simulations in the same manner as used in SDS II. 
We find no clear sign of inhomogeneity in either of (l)-(3) 
for the SDF data. On the other hand, we find a systematic 
difference in the detection limits for the SXDF data with 
a 0.2 magnitude level on arcminute scales. This is proba- 
bly because the majority of the SXDF data were taken in 
2000, when the CCDs installed were a mixture of products 
of two companies and so had a large variety in quantum 
efficiency. Figures 2-5 show, however, that there is no cor- 
relation between the source distribution and the CCD po- 
sitions. This demonstrates that the inhomogeneities seen 
in Figures 2-5 are real, although distributions of LBGs in 
the SXDF on arcminute scales could be contaminated by 
the somewhat large inhomogeneities in the limiting mag- 
nitude. 



3.3. Angular Correlation Function 

In order to quantitatively measure the inhomogeneity 
of spatial distribution, we derive the angular two-point 
correlation function u>(6). According to Landy & Szalay 
(1993), the angular two-point correlation function is cal- 
culated by: 

u obs (6) = \DD{6) - 2DR{6) + RR(0)]/RR(0), (6) 

where DD(6), DR(8), and RR(9) are numbers of galaxy- 
galaxy, galaxy-random, and random-random pairs normal- 
ized by the total number of pairs in each of the three sam- 
ples. We create the random sample composed of 100,000 
sources with the same geometrical constraints as of the 
data sample. The formal error 13 in u>(9) is described by 



y/[l + w ob3 (6)]/DD 



(7) 



(Hewctt 1982). Its 2a error is comparable to the error ob- 
tained with a bootstrap resampling of the data (Baugh, 
Gardner, Frenk, & Sharpies 1996). The real correlation 
function u>(6) is offset from the observed function by the 
integral constant (IC: Groth & Peebles 1977) as, 

u{0) = uj obs (6) + IC. (8) 

We apply the correction for the IC. 14 

Figure 6 shows the angular correlation functions, u>(0), 
of BRi-LBGs (top panel), Viz-LBGs (middle panel), and 
i?iz-LBGs (bottom panel). We measure the angular cor- 
relation function less than ~ 15 arcmin scales, since un- 
certainties in the measurements increase largely over ~ 15 
arcmin scales, which are comparable to a half of the image 
size. In Figure 6, the filled circles are for the SDF and the 
open circles in the top panel are for the SXDF. We find sig- 
nificant clustering signals for BRi-LBGs in both the SDF 
and the SXDF. The data points of the SDF and SXDF 
samples agree within the 2a level except for the points at 
9 ~ 30" and ~ 140", where the measurements of the SXDF 
are significantly lower than those of the SDF. The reason 
for these discrepancies is not clear. However, it could be 
related to the worse quality of the SXDF data than of the 
SDF data. In any case, the SDF sample is deeper than 
the SXDF sample and does not suffer from any detectable 
(artificial) inhomogeneity in the data. Hence, we believe 
that the angular correlation measured for the SDF sam- 
ple is more reliable. We also find clear clustering signals 
for I^iz-LBGs in the SDF. There is a marginal signal for 
Riz-LBGs in the SDF, and no signal is detected for Viz- 
LBGs and Riz-LBGs in the SXDF. This is mainly because 
the numbers of objects in these samples are considerably 
smaller than those in the other samples. 
We fit a single power law, 

u{6) = AZ"6-f>, (9) 

to the data points. For BRi-LBGs, a fit is made to a com- 
bination of the data points of the SDF and the SXDF (A fit 
to the SDF data alone gives almost the same results since 
the errors in the data points of the SDF sample are much 
smaller). Errors in A r ™ and (3 are found to be very large 
except for BRi-LBGs, [3 = 0.90±%m and A™ w = 2.3±° 9 6 

13 The formal error does not include a sample variance, which is caused from field-to-field variations. 

14 If we assume (3 = 0.8, IC/A T J™ is estimated to be ~ 0.006 for our LBG samples, where A[f w and ft are the correlation amplitude and the 
power-law index shown in eq. (9). 
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(Table 3, line 1). This is probably due to the large sam- 
ple size of BRi-LBGs. The value for (3 for BRi-LBGs is 
close to the 'fiducial' value, 0.8. Thus, we use (3 = 0.8 in 
the rest of this paper. Adopting (3 = 0.91 changes the re- 
sults very little. Lines 2-14 of Table 3 present the best fit 
values of A™ w with (3 = 0.8. The correlation amplitudes 
A™ w of BRi-LBGs and Viz-LBGs are as high as at the 
~ 10(7 and ~ 2.5(7 levels, respectively. Hence, their clus- 
tering signals are significant. On the other hand, the A r ™ 
value of Ri z-LBGs is only at the < la significance level. 
Thus, in the following discussion, we adopt the A™ w value 
of SDF Viz-LBGs for the clustering amplitude of z = 5 
LBGs. 

Foreground contamination to a galaxy sample dilutes 
the apparent clustering amplitude of galaxies. When the 
fraction of contaminants is / c , the apparent A™ w value 
can be reduced by a factor of up to (1 — / c ) 2 . The true 
correlation amplitude, A w , is given by 

^raw 

A " = JT^W (10) 

This is the maximum reduction of the correlation ampli- 
tude that occurs when the contaminants are not at all clus- 
tered. In reality, the contaminants in our sample will be 
the sum of foreground galaxies at various rcdshifts, and 
thus would be clustered very weakly, if any, on the sky. 
Thus we use eq. (10) to compute A u for our LBG sam- 
ples. We use f c obtained by simulations described in SDS 
V for LBGs. Table 3 also gives A u values. 

The effect of field-to-field variations in our samples is 
probably modest for LBGs, since we find that angular cor- 
relation functions of BRi-LBGs in the SDF and the SXDF 
agree moderately well within error bars (Figure 6). This is 
because our LBG samples probe large comoving volumes; 
1.6 x 10 6 /if 3 Mpc 3 for BRi-LBGs and 1.7 x 10 6 /if 3 Mpc 3 
for Viz-LBGs. On the other hand, the surveyed volume of 
LAEs at z = 4.9 shown in SDS II is only 9.0 x 10 4 /if 3 Mpc 3 . 
We may have to consider the cosmic variance for the LAE 
sample. 

Since we have a large number of LBGs at z = 4 in the 
SDF, we make subsamples which are divided by magnitude 
or color. We use two observational quantities, i' magnitude 
and i' — z' color, and make three types of subsamples; a 
subsample composed of objects with M < M* , subsamples 
divided by M, and subsamples divided by E(B — V). Then 
we calculate the angular correlation functions of these sub- 
samples in the same manner as for the whole LBG sam- 
ples shown above. We describe the details of the analysis 
in subsections 3.5 and 3.6, and summarize the selection 
criteria and the results in Table 3. We do not apply this 
analysis to the SXDF BRi-LBG sample, since its limiting 
magnitude is brighter than the SDF sample's (and so the 
size of the SXDF sample is about half that of the SDF 
sample) and the SXDF sample could suffer from the inho- 
mogeneity of the limiting magnitude over the image. 

3.4. Correlation Length 

Since the angular correlation function shows clustering 
properties of galaxies projected on the sky, u>(0) reflects 
a combination of the rcdshift distribution of the selected 



galaxies and the intrinsic clustering in the three dimen- 
sional space, i.e., the spatial correlation function, £ s . The 
spatial correlation function of galaxies is usually expressed 
by a power law as 



C 9 = (r/ro)" 



(11) 



where r is the spatial separation between two objects, ro is 
the correlation length, and 7 is the slope of the power law. 
The correlation length, ro, is related to the correlation 
amplitude, A u , with the integral equation called Limber 
equation (Peebles 1980; Efstathiou et al. 1991), 



= Crl / 
Jo 



F{z)Dl^{z)N(zfg(z)dz 



f 

Jo 



N(z)dz 



(12) 

where F(z) 15 describes the redshift dependence of £(r), 
Dg(z) is the angular diameter distance, N(z) is the redshift 
distribution of objects, g(z) = H /c[(l + z) 2 (l + Vl m z + 
Oa((1 + z)~ 2 — l)) 1 / 2 ], and C is a numerical constant, 
C = V^r[( 7 - l)/2]/r( 7 /2). The slope (3 of the angular 
correlation function is related to 7 by 



7 = /?+l- 



(13) 



We need the rcdshift distributions of our LBGs to derive 
r . We adopt the ones obtained in section 2. 

The correlation lengths thus obtained are summarized 
in Table 4. In Table 4, we show both ro and rp aw , which 
are the correlation lengths obtained from A u and A r ™. 
In the following discussion, we regard ro as the best es- 
timate value. We find r = 4.1+oJ and 5 Jt^/iT/ooMpc 
for all LBGs at z — 4 and 5, respectively. Ouchi ct al. 
(2001) calculate the correlation length of LBGs at z = 4 
in the same manner as above, but using a sample con- 
structed based on different selection criteria. They report 
that LBGs at z = 4 with i' < 26 have r = S^t^/i^Mpc 
(r£ aw = 2.7±q j^Mpc). This value is consistent with 
that obtained from our LBG sample. The correlation 
lengths we obtain (r = 3 — 8/iT/qq Mpc) correspond to pro- 
jected angular scales of 100" - 300" at z = 4 - 5. These 
scales fall in the range where the power-law fitting is ex- 
cellent (Figure 6). 

3.5. Evolution of Correlation Length and Galaxy 
Distribution Bias 

We discuss the redshift evolution of the correlation 
length. Since the clustering amplitude (i.e., correlation 
length) depends on the brightness of galaxies (see section 
3.6.1), we compare the correlation lengths of LBGs whose 
luminosities are L > L* . We calculate the correlation 
length of LBGs with L > L* for z = 4 and 5. We make a 
subsample composed of 357 LBGs at z = 4 with i' < 25.3 
corresponding to M < —20.8 (~ M*), and calculate the 
correlation length to be ro = 5.1T^'5^i~oo Mpc (line 2 of 
Table 5). Since the limiting magnitude of z = 5 LBGs is 
z' = 25.8 which corresponds to M = -20.5 (~ M*), we 
regard the results of all z = 5 LBGs (line 2 of Table 4) as 
those of L > L* LBGs (line 3 of Table 5). Thus the corre- 
lation length of L > L* LBGs at z — 5 is r — 



15 Assuming that the clustering pattern is fixed in comoving coordinates in the redshift range of our sample, we take the functional form, 
F(z) = (1 + z)/(l + z c )~(' i+e \ where z c is the central redshift of the sample and e = —1.2. The effect of the change in e over < e < — 3 on 
ro is, however, very small. 
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Mpc. We also show in Table 5 the correlation lengths of 
LBGs with L > L* at z = 3 given by Giavalisco & Dick- 
inson (2001) and of all LAEs at z = 4.9 obtained by SDS 
II. 

The top panel of Figure 7 plots r as a function of red- 
shift, together with r of the underlying dark matter calcu- 
lated by the non-linear model of Peacock & Dodds (1996). 
The correlation lengths of LBGs with L > L* is almost 
constant around r ~ 5h^Q Mpc at z = 3 — 5, although 
it might increase slightly with redshift. It is also found 
that the correlation amplitudes of LBGs (and LAEs) are 
much higher than that of the underlying dark matter. We 
define the galaxy-dark matter clustering bias of the large 
scale (= 8/ifooMpc) as 

= / £ g (r = SforooMpcT" 

9 V&M(r = 8fcra,Mpc)' 

where £,dm is the two-point correlation function of the un- 
derlying dark matter. 

We also calculate r of the dark matter using linear the- 
ory (e.g., Peacock & Dodds 1994), which is plotted in Fig- 
ure 7 by a dashed line, and find a negligible difference from 
the calculation based on the non-linear model. This is be- 
cause the effect of non-linearity is very small on scales of 
r > ro. Jenkins et al. (1998) show that linear theory traces 
the results of TV-body simulations above ro at z = where 
the non-linear effect is larger than at high redshifts. 

Using linear theory, we calculate b g by eq. (14) from ro 
(cq. 11), 7 (eq. 13), and £dm- We explicitly present the 
relation between r and b g given in cq. (14): 

= / [(S/HooMpcVroFT 
9 V&M(r = 8fcra,Mpc)' 

We show values of b g in Table 5 for LBGs with L > L* 
and all LAEs (and Table 4 for all samples), and the bot- 
tom panel of Figure 7 plots b g against redshift for LBGs 
with L > L* and all LAEs. These LBGs and LAEs are 
biased against dark matter by b g <~ 3 — 5, and the bias be- 
comes stronger at higher redshifts. This increase in LBGs' 
bias can be regarded as a piece of evidence supporting the 
biased-galaxy formation scenario (e.g., Baugh et al. 1999). 

3.6. Segregation of Clustering Amplitudes 

3.6.1. Clustering Segregation with UV- Continuum 
Luminosity 

Since we have a large number of LBGs at z = 4, we ex- 
amine the luminosity dependence of clustering amplitude 
(ro and b g ) using four subsamplcs binned according to the 
source brightness (Am = 1.5). We calculate the angular 
correlation functions for the four subsamples, and estimate 
ro and b g in the same manner as described in sections 3.4 
and 3.5, but applying to each subsample its own redshift 
distribution which is determined by simulations described 
in SDS V (Figure 12 of SDS V). Lines 6-9 of Table 3 and 
lines 4-7 of Table 4 give the best-fit parameters of angular 
correlation function together with ro and b g for these four 
subsamples. Figure 8 plots ro and b g for the four subsam- 
ples as a function of magnitude. It is found from Figure 
8 that brighter LBGs are clustered more strongly, which 

16 The relation between E(B — V) and i' — z' is given by E(B — V) 



is a tendency similar to that found for present-day galax- 
ies (Norberg et al. 2002) and LBGs at z — 3 (Giavalisco 
& Dickinson 2001). It is also found that this luminosity 
dependence of ro is weaker for fainter galaxies. We fit the 
data points by a linear function of luminosity (L): 

b g /b* g = a+(l-a)L/L*, (16) 

where a is a free parameter, b g is the galaxy-dark matter 
bias, and b*, which is also a free parameter, is the bias 
of galaxies with luminosity, L* (Norberg's law; Norberg 
et al. 2002). We obtain b* = 2.8 and a = 0.58 with L* 
being fixed at M-j* 700 = —21.0 (or i* = 25.1), which is ob- 
tained in SDS V. The value of a for z = 4 LBGs is smaller 
than that for present-day galaxies, a — 0.85 (Norberg et 
al. 2002). It may indicate that the clustering amplitude of 
z = 4 LBGs depends more strongly on luminosity, and/or 
that the luminosity dependence of clustering is different 
between rest-frame UV luminosity (for LBGs) and optical 
luminosity (for present-day galaxies). 

3.6.2. Clustering Segregation with Dust Extinction 

Strong clustering has been reported for SCUBA sources 
(r - 12.8 ± 4.5 ± 3.0^ o M Pc; Webb et al. 2003). Since 
SCUBA sources are thought to be dust rich starburst 
galaxies, the clustering amplitude of LBGs as a function 
of dust extinction may give a hint for the relation between 
SCUBA sources and LBGs. We calculate the clustering 
amplitude for subsamples of LBGs at z = 4 binned with 
E[B -V). We estimate E(B - V) of LBGs from their 
colors with equations shown in SDS V. We make these 
subsamples from 357 bright LBGs with i' < 25.3 (~ M*) 
so that we accurately measure colors, i' — z' , and thus ob- 
tain E(B — V). We make five subsamples whose E(B — V) 
values range from —0.2 to 0.6 corresponding to i' — z' color 
from -0.18 to 0.50. 16 Lines 10-14 of Table 3 and lines 
8-12 of Table 4 present the best-fit parameters of angular 
correlation functions, together with ro and b g , for these 
five subsamples. We find that the most dusty subsample 
(0.2 < E(B - V) < 0.6, i.e., E(B - V) ~ 0.4) has the 
strongest clustering: ro = 8.2+ 4 '2/ij~ o Mpc. Although the 
errors are large, it may suggest that very dusty LBGs have 
a connection to SCUBA sources. 

On the other hand, Ouchi, Yamada, Kawai, & Ohta 
(1999) find that no LBG identified in the HDF-N is de- 
tected by deep SCUBA observations (Hughes et al. 1998). 
There are three possible explanations for this finding. 
First, dusty LBGs bright enough to be detected with 
SCUBA may be so rare that the probability of finding 
such an LBG in a small area like the HDF-N is very 
low. Second, the redshift distributions of LBGs and 
SCUBA sources may overlap little with each other. In- 
deed, SCUBA sources are generally located at z = 1.9 — 2.8 
(Chapman, Blain, Ivison, & Smail 2003), while LBGs are 
found at redshifts higher than z ~ 2.5. Third, LBGs 
may have a closer relation with EROs than with SCUBA 
sources. EROs are thought to be a mixture of old galax- 
ies and dusty star- forming galaxies. Although we assume 
that the red i' — z' color of LBGs is due to dust extinction, 
the reddest LBG subsample might include a large number 
of old galaxies which have optical-to-infrared colors close 
to those of EROs. In reality, McCarthy et al. (2001) find 

~ 0.0162 + 1.18(i' - z') (see SDS V). 
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functions (LFs) given in SDS V. Note that, since the LFs 
are derived by correcting completeness and contamination 
for the number counts of LBGs, these number densities do 
not suffer from these systematic observational errors. The 
ticks in Figure 8 show the halo masses and the number 
densities for z = 4 LBGs thus obtained. We then plot in 
Figure 9 the number density of LBGs against the mass of 
hosting dark halos for each luminosity bin, together with 
the mass function of all dark halos existing at z — 4 cal- 
culated by eq. (17). 

4.2. Hosting Dark Halos and Halo Occupation Number 

4.2.1. Dark Halos of Lyman Break Galaxies at z = 4 

Figure 9 shows that dark halos which host LBGs at 
z = 4 have 10 11 - 10 13 h^M Q . Each filled circle denotes 
the mass of hosting halos for each luminosity bin (see Ta- 
ble 6 for the definition of the luminosity bins); from right 
to left, halos hosting LBGs with i' ~ 24.0, 24.5, 25.0, 
and 25.5. The third circle from right is for i' ~ 25.0 
(Mi 700 ^ -21) LBGs, i.e., for L*-LBGs (sec SDS V). 17 
We find that bright LBGs (M ~ M* — 1) have more mas- 
sive dark halos {4.Q±f% x W 12 h^M Q ) than faint LBGs 

(M ~ M* + 0.5; 1.7to.7 x lO 11 /^ 1 -^©)- A qualitatively 
similar trend has been found for LBGs at z = 3 (Giavalisco 
& Dickinson 2001). 

The most striking feature in Figure 9 is the discrepancy 
in number density between all dark halos (thin solid line) 
and all the data points (filled circles) except the brightest 
one. We define the occupation number, -/V occup , as the ra- 
tio of the number density of galaxies to that of dark halos. 
Table 6 presents iV occup for each luminosity bin together 
with the results for galaxies at z — 3 and 5, which are 
explained in the following sections. Lines 10-13 of Table 
6 (and Figure 9) show that the occupation number of the 
brightest LBGs is almost unity, but those of fainter LBGs 
(including L*-LBGs) are 3 to 10 times smaller. This im- 
plies that the majority of low-mass (< 10 12 h? Q Mq) dark 
halos have no LBGs. Since LBGs are star-forming galaxies 
which are bright in the UV, it is possible that low-mass 
dark halos have a galaxy which happens to be between 
star-forming phases and thus escapes from our LBG sam- 
ples. 

4.2.2. Dark Halos of Galaxies at z = 3 

In order to investigate the discrepancy in number den- 
sity between LBGs and all dark halos (Figure 9), we 
carry out the same analyses for various galaxy populations 
at z ~ 3, whose number density and correlation length 
have been reported in the literature. We examine K- 
band selected galaxies (FIRES; Daddi et al. 2003), SCUBA 
sources (SCUBA; Webb et al. 2003), and LBGs (Giavalisco 
& Dickinson 2001). Here, we note basic properties of 
these three populations. FIRES galaxies are detected by 
very deep X-band imaging (K < 24) with VLT/ISSAC 
(Daddi et al. 2003). Although the sample of FIRES galax- 
ies may suffer from the field variance due to its very small 
survey field (~ 4 arcmin 2 ), Daddi et al. (2003) claim from 
their simulations that the strong correlation amplitudes 

17 We also have a subsample of L > L* LBGs which contains all LBGs with L > L* (see section 3.3). Note that the samples of L*-LBGs and 
L > L* LBGs arc different samples. 



that the correlation length of EROs at z <~ 1 is 9 — 10h 1QQ 
Mpc, comparable to that of the reddest LBG subsample. 

4. GALAXY AND STRUCTURE FORMATION BASED ON THE 
HIERARCHICAL CLUSTERING PICTURE 

4.1. Joint Analysis of Luminosity Functions and 
Clustering Amplitudes based on the Cold Dark 
Matter Model 

We estimate the mass of dark halos which host LBGs 
and LAEs using the analytic models given by Sheth & 
Tormen (1999) and Sheth, Mo, & Tormen (2001). They 
are derived from a fit to the results of GIF simulations 
(GIF/ Virgo collaboration; e.g., Kauffmann, Colberg, Di- 
aferio, & White 1999). 

Prior to estimating the mass of dark halos, we present 
an overview of the model of Sheth & Tormen (1999). Ac- 
cording to their model (see also Mo & White 2002), the 
number density of collapsed dark halos with mass M at 
redshift z is given by 

n(M,z) d M = A(l + ^)/f |^exp d M, 

(17) 

where v' = y/av, a = 0.707, A ~ 0.322, and q = 0.3. Here, 
v is defined by 

V ^ D(z)a(M) ' (18) 

where D(z) is the growth factor, a(M) is the rms of the 
density fluctuations on mass scale M, and S c = 1.69 rep- 
resents the critical amplitude of the perturbation for col- 
lapse. We calculate the growth factor following Carroll, 
Press, & Turner (1992). a(M) is calculated from the 
power spectrum whose power law index is n = 1, using 
the transfer function given by Bardeen, Bond, Kaiser, & 
Szalay (1986). The number density of collapsed dark halos 
at a given redshift, i.e., the mass function of dark halos, 
is calculated using eqs. (17) and (18). The resulting mass 
function at z — 4 is shown in Figure 9. 

In the model of Sheth & Tormen (1999), the bias of dark 
halos against the dark matter, 6dh, is calculated by 



3DH 



1 



.'2 



+ bv' 2 ^ - 



/2c 



v' 2c + &(l-c)(l-c/2) 

(19) 

where b = 0.5, c = 0.6, and the other parameters (z/, <5 C , 
and a) are the same as those in eqs. (17) and (18) (Sheth, 
Mo, & Tormen 2001). Since v' is a function of redshift 
and mass, the mass of a dark halo is estimated from eqs. 
(18)and (19), once the bias of the dark halo and the red- 
shift are given. Since the bias values of LBGs, b g , are 
measured at a large separation of 8/ir/ o Mpc, these values 
reflect the bias of dark halos (which host LBGs) rather 
than the bias of the galaxy distribution in dark halos, i.e., 
b g ~ &dh- Assuming that galaxies with a given brightness 
are hosted by dark halos of a specific mass (i.e., there is 
a one-to-one correspondence between the galaxy luminos- 
ity and the halo mass), we estimate the average mass of 
dark halos, (M), for galaxies of a given luminosity range 
from the bias values obtained in section 3.5. We also cal- 
culate the number density of galaxies from the luminosity 
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they found are hardly reproduced by the field variance 
alone. FIRES galaxies include UV- faint galaxies, i.e., qui- 
escent in star-formation activity, which are not identified 
by LBG selections (Franx et al. 2003). SCUBA sources are 
detected by their strong submillimctcr fluxes in Canada- 
UK Deep Submillimeter Survey with JCMT/SCUBA, and 
they are thought to be a mixture of dusty-starburst galax- 
ies and AGNs. The sample of LBGs at z = 3 used in this 
analysis is taken from Giavalisco & Dickinson (2001) (see 
also Giavalisco et al. 1998; Adelberger et al. 1998; Arnouts 
et al. 1999, 2002; Porciani & Giavalisco 2002). 

Figure 10 (and lines 1-7 of Table 6) shows the results 
for z ~ 3 galaxies. The mass of dark halos hosting LBGs 
ranges over 10 11 — lO 12 /^ M Q . The occupation number 
of LBGs is less than unity (<~ 0.3: see Table 6). These 
features of LBGs at z = 3 are similar to those at z = 4. 
On the other hand, dark halos of FIRES galaxies have 
~ lO 12 /i7Q 1 M , which is comparable to the mass of LBGs, 
but the occupation number of FIRES galaxies is as high 
as ~ 4. This comparable mass and large occupation num- 
ber indicate that dark halos with ~ 10 12 ft-f 1 M Q have 4 
galaxies on average, and only about 1/10 (~ 0.3/4) of 
them have an active-star formation whose star-formation 
rate (SFR) exceeds ~ S/i^Mgyr" 1 (or ~ 20/if 2 M Q yr" 1 
if dust extinction is corrected; see SDS V), so that they 
arc identified as LBGs by their bright UV continuum. 

Halos hosting SCUBA sources are located near the high- 
mass end of the dark halo mass function at z = 3, having 
~ 3 x lO 13 /).^ M Q . The best estimate of the occupation 
number of SCUBA sources is about 4, though the uncer- 
tainty in it is about an order of magnitude larger than 
those for other objects due to large errors in the measure- 
ments of the number density and the correlation length. 
The large mass of dark halos (~ 5 x l0 13 hjQ Mq) hosting 
SCUBA sources indicates that SCUBA sources exist only 
in proto-cluster regions. 

4.2.3. Dark Halos of Lyman Break Galaxies and Lyman 
a Emitters at z = 5 

We then apply our analysis to LBGs at z = 5 and LAEs 
at z = 5 obtained in SDS II, so as to investigate the evolu- 
tion of LBGs from z = 3 to z = 5 and differences between 
LBGs and LAEs. The results are given in lines 14-15 of 
Table 6 and in Figure 11. Both LBGs and LAEs at z = 5 
are embedded in dark halos with <~ 10 12 H^qMq. However, 
the number densities of LBGs and LAEs are different, and 
so the occupation number of LAEs (N occup = 3.2) is about 
5 times larger than that of LBGs (A occup = 0.6). It indi- 
cates that the majority of dark halos with ~ lO 12 /!^ M@ 
do not host LBGs, while hosting multiple LAEs. Since the 
typical brightness of LAEs in SDS II is about 2 magnitude 
fainter in UV continuum than that of our LBGs at z = 5, 
galaxies in the LAE sample have about ten times less star- 
formation rates than those in the LBG sample. This deep 
detection limit of the LAEs results in the larger number 
density for LAEs. 

Note that the estimated mass of dark halos hosting 
LAEs may be smaller than the true values. We use 
the best-fit power-law function for the angular correla- 
tion function shown in SDS II to obtain r , and calculate 
the galaxy-dark matter bias b g defined at r — 8/iJ~ 00 Mpc. 
However, the angular correlation function has a hump at 



around r = 8h^ Mpc, and the observed points systemati- 
cally exceed the fitted power law (see Figure 6 of SDS II). 
If the hump of the clustering amplitude at r ~ 8/iJ~ 00 Mpc 
is real, the true b g value is larger than that obtained 
in section 3.4, resulting in a larger mass of dark halos 
(3.5^07 x 10 12 IijqMq) and a larger occupation number 
(35). ■ 

4.3. Descendants of Dark Halos Hosting High-z Galaxies 

Once the mass of dark halos is determined at a given 
rcdshift (z), one can calculate analytically the bias of dark 
halos for their descendants (and progenitors) at an ar- 
bitrary redshift. The bias of descendants at the present 
epoch, &q H , is calculated by (Sheth, Mo, & Tormen 2001): 



1+ 



D(z) 



v n + &! /2(i- c ) 



,11c 



+ &(l-c)(l 



c/2). 
(20) 

where D(z) is the growth factor, and the other notations 
are the same as those in eq. (19). Assuming that the bias 
thus obtained represents the bias of a given dark-halo mass 
at the present epoch, we estimate the dark-halo mass of 
descendants with eq. (19). We first calculate the present- 
day bias for dark halos hosting high-z galaxies, and then 
derive the mass of descendants at z = from 6p H . 

Figure 12 shows the mass of descendants of dark halos 
which hosted high-z galaxies. In this figure, the num- 
ber densities of descendants are the same as those mea- 
sured for high-z galaxies, since the analytic model does 
not predict the number density of descendants. In other 
words, in Figure 12 we assume that high-z galaxies con- 
serve the number densities until the present epoch. It is 
found from Figure 12 that the majority of high-z galax- 
ies are embedded in dark halos of 10 13 — 10 15 /i 70 1 M Q at 
the present epoch. Generally speaking, dark-halo masses 
of 10 13 — 10 15 /i 70 1 Mq at z = correspond to the mass 
scale of clusters and groups. Our finding means that the 
high-z galaxies examined here become member galaxies of 
clusters and groups in the present-day Universe in a sta- 
tistical sense. This result is qualitatively consistent with 
those found by Nagamine (2002) who claims from a hy- 
drodynamic simulation that all bright (My < —23) LBGs 
at z = 3 evolve into galaxies in clusters or groups at z = 0, 
and that half of LBGs with -23 < M v < -21 fall into 
clusters or groups. 

In Figure 12, the dotted line denotes the number den- 
sity of present-day galaxies, which is calculated from the 
occupation-number function derived by van den Bosch, 
Yang, k Mo (2003) who used the 2dFGRS data. Note that 
this occupation-number function is for galaxies brighter 
than Mb ~ — 19 (this limiting magnitude is determined by 
the depth of the 2dFGRS data). In other words, the dot- 
ted line of Figure 12 indicates the number density of nor- 
mal bright galaxies. This number density of present-day 
normal galaxies is found to exceed those of bright LBGs 
and SCUBA sources. On the other hand, the number den- 
sity of faint LBGs (L < L*), FIRES galaxies, and LAEs 
is higher than that of present-day normal galaxies. If we 
assume that all high-z galaxies evolve into present-day nor- 
mal galaxies, the excess of the number density of high-z 
galaxies over that of present-day galaxies puts a constraint 
on the minimum number of mergers which high-z galaxies 
have experienced until z = 0. We define 7V mm (merge) as 
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the ratio of the number density of high-z galaxies to that of 
present-day galaxies, and calculate 7V min (merge) for each 
of the high-z galaxy populations in Table 6. It is found 
that L < L* LBGs, LAEs, and FIRES galaxies should ex- 
perience mergers at least 2-5, 4, and 7 times, respectively, 
so as to reduce their number density. Part of these high-z 
galaxies will evolve into massive cluster galaxies through 
mergers. Our findings are consistent with results of semi- 
analytical models. For example, Governato et al. (2001) 
find that among 12 halos containing at least one LBG in a 
proto-cluster region at z = 3, 7 halos merge into one cen- 
tral object of the cluster while the rest of 5 halos survive 
as separate entities inside the cluster at the present epoch. 

The thick solid curve in Figure 12 presents the relation 
between the observed number density of z = 4 LBGs and 
the predicted dark-halo mass of their descendants at z = 0. 
This curve is derived by evolving the thick solid curve in 
Figure 9 to the present epoch using eq. (20). It is interest- 
ing to note that the data points for LBGs at z = 3 and 5 
are located close to this solid curve. This may imply that 
descendants of LBGs at z = 3, 4, and 5 are indistinguish- 
able in terms of the number density and halo mass. 

We find that the mass of descendants of dark halos host- 
ing SCUBA sources is about 4.51^ x IO^^Mq. Thus 
the descendants of SCUBA sources will probably be em- 
bedded in clusters of galaxies. We then find that the data 
point of SCUBA sources is located near the high-mass 
end of the thick solid curve in Figure 12. This would 
indicate that SCUBA sources are very dusty and high 
star-formation rate LBGs (a similar discussions is found 
in Adelberger & Steidel 2000). The number density of 
SCUBA sources is about 1/10 that of present-day normal 
galaxies in clusters. This implies that the descendants of 
SCUBA sources are a rare population in present-day clus- 
ter members. Since the fraction of elliptical galaxies in 
present-day cluster members is about 20% (Dressier 1980; 
Whitmore, Gilmore, & Jones 1993), the descendants of 
SCUBA sources could be cluster ellipticals. This infer- 
ence supports the prediction by Shu, Mao, & Mo (2001) 
based on the CDM model that descendants of bright sub- 
mm sources should reside in present-day clusters and that 
these objects are likely to be the progenitors of giant el- 
lipticals. 

5. DISCUSSION 

5.1. Clustering Properties of LBGs at z — 4 and 5 

Figure 7 shows the observed correlation length, r (in 
comoving units), and galaxy-dark matter bias, b g , of LBGs 
at z ~ 4 and 5 in our sample, together with those of galax- 
ies at various redshifts between z — and 3 taken from 
the literature. The data points for LBGs (at z = 3 — 5) 
are for galaxies with L > L* in UV continuum, while the 
other points at < z < 1 are for galaxies selected by op- 
tical continuum. Figure 7 indicates that the clustering of 
LBGs at z = 3 — 5 is strongly biased against the under- 
lying dark matter (dotted or dashed line). The values for 
the bias parameter are b g ~ 3 — 5 at z = 3 — 5. The bias 
increases monotonically with redshift. The increase in b g 
with redshift has been predicted by semi-analytic mod- 
els (Baugh et al. 1999; Kauffmann, Colberg, Diaferio, & 
White 1999) and hydrodynamical simulations (Blanton et 
al. 2000; Yoshikawa, Taruya, Jing, & Suto 2001). The ob- 



served b g values are consistent with predictions by Kauff- 
mann, Colberg, Diaferio, & White (1999) (solid line of Fig- 
ure 7) for galaxies with M < —19 + 5 log h in the B band. 
Although the selection criteria of galaxies arc different be- 
tween ours (LBG selection) and Kauffmann, Colberg, Di- 
aferio, & White's (B-band selection), their semi-analytic 
model reproduces the observed trend of b g fairly well. 

In Figure 7, ro decreases with redshift over z ~ — 1, 
has a minimum at z <~ 1 — 3, and then slightly in- 
creases at z > 3. This behavior can be explained by 
a combination of (i) the increase in bias (b g ) with red- 
shift and (ii) the decrease in the correlation amplitude of 
dark matter ((dm) with redshift, since r is expressed as 
r /(8Mpc) = 4 2/7) fr>M(r = 8Mpc)(V7) ( eq . 15). At low 
redshifts, the decrease in (dm (with redshift) dominates 
over the increase in b g . This is the reason why ro de- 
creases with redshift over z <~ — 1. On the other hand, 
the increase in b g is more effective at high redshifts, result- 
ing in an increase in ro with redshift at z > 3. Note that 
these two factors roughly cancel with each other, and thus 
the change in ro over z = — 5 is modest compared with 
the change in b g . Kauffmann, Colberg, Diaferio, & White 
(1999) have found a qualitatively similar trend in the evo- 
lution of ro in their galaxy formation models. They have 
claimed that a dip in ro between z = and 1 occurs in 
A-dominated flat universes, because structures form early 
and dark halos hosting L* galaxies are unbiased tracers of 
the dark matter over this redshift range. To summarize, 
the observed behaviors of b g and ro found in this paper can 
be regarded as a piece of evidence supporting the biased 
galaxy formation scenario. 

5.2. Dark-Matter Properties of High Redshift Galaxies 

It is found from Figures 9 and 10 that the mass of 
dark halos hosting LBGs at z = 4 — 5 ranges from 10 11 
to 10 13 /if 1 M Q . Giavalisco & Dickinson (2001), Bul- 
lock, Wechsler, & Somerville (2002), and Moustakas & 
Somerville (2002) have estimated the mass of dark ha- 
los hosting LBGs at z = 3 using analytic approaches 
based on the CDM model. Giavalisco & Dickinson (2001) 
have found that the average mass of dark halos for LBGs 
with R = 23, 25.5, and 27 is (M) = 3.6, 1.3, and 
0.6 x 10 12 h? Mq, respectively. Bullock, Wechsler, & 
Somerville (2002) have derived the minimum mass of dark 
halos, M m i n , and the mass of dark halos hosting one 
LBG, Mi, using the observed spatial correlation functions 
given by Adelberger et al. (1998), Giavalisco & Dickin- 
son (2001), and Arnouts et al. (1999). They have found 
M min ~ (0.6 - 1.1) x 10 w h^M Q and M 1 = (0.9 - 1.4) x 
lO 13 ftf o 1 M . Moustakas & Somerville (2002) have ob- 
tained M min ~ 1.9 x 10 w h^M Q , Mi ~ 8.6 x 10 12 h^M e , 
and (M) ~ 8 x lO 11 ft^ 1 M . ^From Figure 10, dark halos 
hosting LBGs at z = 3 have Mi ~ 10 13 h^M Q , which 
is close to the values obtained by Bullock, Wechsler, & 
Somerville (2002). Similarly, the typical mass of dark ha- 
los hosting LBGs at z = 3, (M) ~ 10 12 h^M Q , is close to 
those derived by Giavalisco & Dickinson (2001) and Mous- 
takas & Somerville (2002). Thus, our results for LBGs at 
z = 3 (M = 1 - 10 x 10 12 ) are consistent with those 
given by the previous works. Furthermore, Hamana et 
al. (2004) have derived recently the mass of dark halos 
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of our LBGs and LAEs at z = 4 and 5 using the halo 
occupation function. They directly fit predicted angular 
correlation functions and number density to data for LBGs 
and LAEs (obtained in this paper and SDS II), and obtain 
(M), Mi, and M m ; n . They assume the power-law occupa- 
tion number function and the same analytic halo models 
(Shcth & Tormen 1999) as used in this paper. They obtain 
(M) = 9 x 10 11 , 6 x 10 11 , and 1 x 10 12 h^M Q for LBGs 
at z — 4 and 5 and LAEs at z = 5, respectively, which are 
comparable to our results (M = 4 x 10 11 , 9 x 10 11 , and 
1 x 10 12 ). Hamana et al. (2004) have found an indication 
that the simple halo model may not work for small-scale 
clustering (0 < 120" corresponding to r < 3.1/i7/ o Mpc) 
of LAEs. In our analysis, we obtain the bias at a large 
scale (at 8h^ Mpc), and do not examine whether £dh is 
fit to £(LAE) on small scales. In this sense, our analysis 
is simpler than theirs. 

In Table 7, we summarize the properties of dark halo 
hosting L > L* LBGs (and LAEs), which are extracted 
from Table 6. We find the typical mass of dark halos 
hosting L > L* LBGs is 1.3^ x 10 12 h^M Q , 1.0±° 9 7 x 
10 12 h^M Q , and x 10 11 /if 1 M Q for z = 3, 4, and 

5, respectively. Thus the dark halo masses hosting L > L* 
LBGs are almost constant around 1 x 10 12 h 7 Q M which 
is comparable to the total mass of the Milky Way Galaxy 
at the present epoch. Although there might be a trend 
that the typical mass decreases slightly with redshift, we 
find no clear difference in the mass range of dark halos 
hosting LBGs with L > L* (M < -20.5 or, equivalently, 
SFR> 10 [4O]/i7 o 2 M yr _1 before [after] extinction correc- 
tion) over z = 3 — 5. This finding suggests that the mass- 
to-luminosity ratio of LBGs (the mass of a hosting halo 
divided by a UV-luminosity of the LBG in the halo) does 
not change largely from z — 3 to 5. In other words, the 
star formation efficiency is almost constant over z = 3 — 5. 
This suggestion may not be consistent with the theoretical 
prediction that the star-formation efficiency increases with 
redshift (e.g., Hernquist & Springcl 2003). 

Figure 9 and Table 6 show that the occupation number, 
Aoccup, of LBGs for massive dark halos (<~ 10 13 /i7 M ) 
is close to unity while A occup for less massive halos (~ 
10 12 hyQ M & ) is about 0.1. Thus, it is concluded that more 
massive halos possess a larger number of LBGs. Note that 
a similar feature is found for the 2dFGRS galaxies at z = 
in the range of dark halo mass of of 10 12 — 10 13 M & as 
shown in Figure 12. 

5.3. Unified View for Various Kinds of High-z Galaxies 

In this section, we discuss relations between LBGs, 
LAEs, X-band selected galaxies (FIRES), and SCUBA 
sources. These galaxies are named after the search tech- 
niques (section 1), and are all located at high redshifts. 
Little has been known to date about relations between 
these galaxy populations. 

Figure 12 plots the observed number densities of these 
high-z galaxies against the predicted masses of dark halos 
of their descendants at z = 0. In this figure, the thick 
solid curve corresponds to the best-fit function to the ob- 
served relation between the bias parameter and the num- 
ber density for our z = 4 LBGs (Figure 8 of section 3.6.1). 
Interestingly, the data points of LBGs at both z = 3 and 
z = 5 (blue and red circles) are located very close to this 



curve. This suggests that LBGs at z = 3 — 5 are a single 
population in terms of the mass of hosting halos and the 
halo occupation number. 

We propose here that LBGs at z = 3 — 5 come from 
a single parent population of galaxies, assuming that any 
dark halo more massive than a critical value has such par- 
ent galaxies and that the stellar mass of these galaxies is 
proportional to the mass of the dark halo. These parent 
galaxies can be detected by NIR surveys such as FIRES, 
because NIR surveys sample galaxies according to their 
rest-frame optical luminosity, i.e., their stellar mass. Fig- 
ure 10 shows that the number density of FIRES galaxies 
is 10 times higher than that of LBGs at the same mass 
range of dark halos. If LBGs are star- forming galaxies 
with intermittent star formation activities as suggested 
in section 4.2.1, this difference in the number density re- 
flects probably the total duration time of active star for- 
mation in galaxies. Both LBGs and FIRES galaxies are 
seen at z = 2.5 and 3.5, and the cosmic time between 
these two redshifts is 8.0 x 10 8 yrs. If we assume that 
LBGs are visible in optical observations (i.e., bright in the 
rest-frame UV wavelength) only when they are in a phase 
of active star formation, then the total duration time of 
star formation between z — 2.5 and 3.5 is estimated to be 
8.0 x 10 8 /10 = 8 x 10 7 yrs. This value is comparable to 
the typical age of stellar populations of LBGs, 7 x 10 7 yrs, 
given by Papovich, Dickinson, & Ferguson (2001). They 
have derived the age from a fitting of model spectra to 
the observed optical-to-NIR colors of LBGs at z ^ 3 , as- 
suming that the star formation history of LBGs is ap- 
proximated by an exponentially decaying star-formation 
rate. Furthermore, if we assume that the typical star- 
formation rate of LBGs is 20— 2OOM /yr (whose values are 
extinction-corrected by a factor of 4: see Figure 16 and sec- 
tion 5 of SDS V), then the stellar mass accumulating from 
z = 5 to 3 (At = 9.6 x 10 8 yrs) is estimated to be roughly 
9.6x 10 8 [yr]/10x (2O-2OO)[M yr- 1 ] - (2 - 20) x 1O 9 M . 
This value is comparable to the typical stellar masses of 
LBGs at z ~ 3 derived by Shapley et al. (2001) and Pa- 
povich, Dickinson, & Ferguson (2001). On the other hand, 
if LBGs are FIRES galaxies in a phase of star formation, 
LBGs should have similar stellar masses to FIRES galax- 
ies. This prediction is supported by Franx et al. (2003) 's 
finding that while FIRES galaxies are slightly less lumi- 
nous in the observed-frame K band than the brightest 
LBGs, their stellar masses are comparable to, or higher 
than, the brightest LBGs'. 

The data point of SCUBA sources (star mark) is lo- 
cated at the high-mass end of the solid curve in Figure 
12. This coincidence would indicate that SCUBA sources 
and LBGs belong to the same population. As discussed in 
section 3.6.1, brighter LBGs have larger dark-halo masses, 
indicating that LBGs with higher star-formation rates re- 
side in more massive dark halos. The star formation rates 
of SCUBA sources, ~ 10 3 M Q /yr, are (much) higher than 
those of LBGs, 20 — 2OOM /yr after extinction correction. 
Thus, the correlation found for our LBGs at z = 4 holds 
for SCUBA sources as well. We have also found a possi- 
ble connection of dusty LBGs to SCUBA sources in the 
analysis of section 3.6.2. 

We have found a probable connection between LBGs, 
FIRES galaxies, and SCUBA sources as discussed above. 
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However, LAEs at z = 4.9 cannot be understood in a 
similar way. Since LAEs are also star-forming galaxies, 
they should be located near the solid curve in Figure 12, 
if they belong to the same population as LBGs. In this 
case, LAEs are expected to be located near the low-mass 
end of the solid curve, because the typical star-formation 
rate of LAEs is only about 10 (3) M /yr with (without) 
dust extinction correction. 18 We think that this appar- 
ent discrepancy is related to (i) a large cosmic variance 
in the LAE sample of SDS II and/or (ii) a selection bias 
of LAEs. First, as we note in section 3.3, the surveyed 
volume for their LAEs is just 9.0 x 10 4 /if 3 Mpc 3 , and thus 
the results could be strongly affected by the field variance. 
In reality, Shimasaku et al. (2003) have found that their 
LAEs make a large-scale structure, and it is possible that 
the clustering derived in SDS II is biased high due to the 
large-scale structure. Second, star-forming galaxies with 
relatively bright Lyman a emissions can be selected as 
LAEs even when their star formation rates (i.e., UV lu- 
minosities) are considerably low. We have found in SDS 
II that LAEs brighter in Lyman a emission are clustered 
more strongly while clustering strength does not clearly 
correlate with UV luminosity. If this correlation holds for 
all star-forming galaxies, the clustering strength of LAEs 
should be higher than that of general star-forming galaxies 
with similar star-formation rates, resulting in very massive 
dark halos for LAEs. The cause of this 'Lyman a lumi- 
nosity bias' is not clear to us, but a possible explanation is 
that the Lyman a luminosities of galaxies increase (only) 
when they infall toward massive dark halos. When galax- 
ies infall toward massive halos, star formation may be trig- 
gered in them due to, for instance, interactions with the 
halos. Since the duration of strong Lyman a emission is 
very short, galaxies with strong Lyman a emission will 
be found only around massive halos, and thus LAEs have 
strong clustering. Scannapieco & Thacker (2003) have dis- 
cussed a similar scenario for high-z starburst galaxies. In 
this explanation, we assume that Lyman a luminosity cor- 
relates with star formation rate fairly well. However, this 
assumption may not necessarily be true. For instance, 
Shapley et al. (2001) find that LBGs with stronger Lyman 
a emissions have older stellar populations, contrary to the 
predictions from simple Hn region models. In any case, 
our knowledge about the nature of LAEs is still poor. It 
is of great importance to accumulate the data of LAEs 
by new observations to infer more accurately their funda- 
mental quantities such as halo mass, stellar mass, and star 
formation rate. 

Finally, we summarize our unified view for high-z galax- 
ies (i.e., LBGs, FIRES galaxies, SCUBA sources, and 
LAEs), which we have proposed in this subsection. At 
z ~ 3 — 5, dark halos with ~ 10 12 /iy M host a few galax- 
ies (section 4.2.2) whose stellar masses are large enough to 
be detected as FIRES galaxies. Most of the galaxies hosted 
by these dark halos are quiescent in star-formation, but 
about one tenth have a moderately strong star-formation 
activity (SFR~ 20 - 2OOM /yr), producing a bright UV 
continuum and thus being identified as LBGs. Hence, 



FIRES galaxies are the parent population for LBGs ac- 
cording to our view. LBGs found at z <~ 3,4, and 5 
are not necessarily on the same evolutionary track, since 
the duration of a single star formation for LBGs is ex- 
pected to be much shorter than the cosmic times be- 
tween z — 3, 4, and 5. More massive dark halos tend 
to host LBGs with higher SFRs, and extremely massive 
dark halos with <~ lO 13 /^ M have a starburst galaxy 
with SFR- 10 3 M o /yr, which is identified as a SCUBA 
source. The relationship between LAEs and the other 
high-z galaxies is not clear. 



5.4. Formation History of Galaxies and Large-Scale 
Structures Based on the Observational Results 

In sections 5.3, we have proposed that LBGs and 
SCUBA sources are high-z star-forming galaxies emerged 
from the same parent population. In our proposal, high- 
z galaxies experience many episodic star-forming phases 
with star-formation rate of 20 — 2OO/i 7o 2 M yr _1 , whose 
characteristic total duration time is ~ 80Myr (per 1 Gyr) . 
The diversity of the properties of star-formation galaxies 
may be due to the variety of star-formation rate. These 
galaxies accumulate stars through episodic star-formation 
activities of 20 — 200h^ M yr _1 , and have stellar masses 
of - 1 x 1O 1O M at z = 3. 

In SDS V, we have found that the luminosity function 
of LBGs does not change from z = 4 to z — 3, while 
the number of bright LBGs (M1700 < —22; extinction- 
corrected SFR is > 100/if 2 M Q yr _1 ) drops toward z = 5. 
If we assume that the star-formation efficiency is constant 
from z = 5 to z = 4, this implies that massive galax- 
ies (with large cold gas reservoirs) are rare at z > 5 
since matter (baryon + dark matter) has not been assem- 
bled by the gravitational instability at z = 5. However, 
the cosmic star-formation rate density is almost constant 
from z ~ 5 to z ~ 1 (SDS V). This is because the ma- 
jor contributors to the cosmic star formation are galaxies 
with smaller star-formation rates ~ I/i^q M yr _1 (SDS 



V). Then, most of the galaxies with active star formation 
of ~ lOO/iyo M yr _1 disappear in the present-day Uni- 
verse (Figure 16 of SDS V). This reduction of galaxies 
with intensive star-forming activities leads to the decrease 
in the star- formation rate density from z ~ 1 to z = 0. 

On the other hand, dark matter is continuously assem- 
bled through the cosmic time (Figure 7). Thus the clus- 
tering amplitude, or correlation length, of the dark matter 
monotonically increases with time. However, the correla- 
tion length of galaxies with L > L* docs not significantly 
change from z = 5 to z = because the galaxy-dark mat- 
ter bias (b g ) decreases with time. Thus, the galaxy clus- 
tering at z = 5 is similar to that at the present epoch. 
This is the reason why SDS IV discovered the large-scale 
structure of galaxies at z = 5. The size of the large-scale 
structure at z = 5 is comparable to those in the present- 
day Universe found by, e.g., Geller & Huchra (1989). 

We have found that dark halos hosting high-z galax- 
ies (LBGs, LAEs etc.) are considerably massive, ~ 

18 We estimate these SFR values from their median UV continuum flux. If we assume LAEs to be dust free, we obtain the typical SFR of 
about 3M0/yr. On the other hand, if we assume that the dust extinction of LAEs is comparable to that of LBGs at z=4, i.e., E(B-V)=0.15, 
we obtain IOMq/jt for the extinction corrected SFR. These values are comparable to, or lower than, the average extinction-corrected SFR for 
the faintest subsample of BRi-LBGs, although the dust extinction amount of LAEs is not yet well constrained. However, as found from Figure 
12, LAEs deviate largely from the solid curve toward higher halo masses; dark halos hosting LAEs arc as massive as those hosting bright LBGs. 
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10 12 Ii^qMq, which is comparable to the total mass of the 
Milky Way Galaxy. The values of the bias parameter for 
such massive halos are very high at the observed rcdshifts. 
The large masses of dark halos hosting high-z galaxies im- 
ply that they will have evolved into groups and clusters of 
galaxies at the present epoch (section 4.3 and Figure 12) in 
a statistical sense. This would indicate that the majority 
of present-day field galaxies were formed at z < 3, or that 
high-z progenitors of field galaxies are much fainter than 
the limiting magnitudes of today's observations (i' ~ 27 
and K ~ 24). 

6. CONCLUSIONS 

We detect 2,600 Lyman break galaxies (LBGs) at 
z = 3.5- 5.2 in the deep (i' ~ 27) and wide-field (1,200 
arcmin 2 ) data in the Subaru Deep Field (SDF) and the 
Subaru/XMM Deep Field (SXDF). First, we derive clus- 
tering properties (angular correlation function and corre- 
lation length etc.) of LBGs at z — 4 and 5 and present 
the observational results. We then combine these ob- 
servational results with the Cold Dark Matter model to 
understand the formation and evolution of galaxies (lumi- 
nous matter) and dark halos (dark matter) simultaneously 
and obtain theoretical implications. The major findings 
of our study are categorized in (a) observational results 
and (b) theoretical implications as follows. 

(a) Observational results: 

1 . We calculate angular two-point correlation functions 
of LBGs, and find clear signals for LBGs at z = 4 and 5. 
We estimate the correlation length, tq, of spatial two-point 
correlation function (£ g ) by Limber deprojection using the 
redshift distribution obtained by simulations (section 2). 
We obtain r = 4.1 ± 0.2h^ Mpc and 5.9t\ 3 h^ Mpc 
for LBGs at z = 4 and z = 5, respectively. We also cal- 
culate r for bright galaxies (L > L* , i.e., i' < 25.3 for 
LBGs at z = 4 and z' < 25.8 for LBGs at z = 5) to ob- 
tain r = 5.11^/1^0 Mpc and 5.9t^/i^ o Mpc for LBGs 
at z — 4 and 5, and compare them with that of z = 3 
LBGs with similar luminosities, to find that the correla- 
tion length is almost constant at <~ 5hi 00 Mpc from z = 3 
to z = 5. 

2. We find that the correlation amplitude of galaxies is 
larger than that of underlying dark-matter ((dm) at high 
rcdshifts (at least at z > 3). We calculate the galaxy-dark 
matter bias, b g = a/£ s /£dm at 8/if oMpc, for LBGs. We 
find that b g of L > L* LBGs at z = 4 and z = 5 is 3.5±g;f 

and 4.6^1.2) implying that the distribution of LBGs at 
z > 3 is highly biased against the underlying dark mat- 
ter. The bias of galaxies brighter than ~ L* monotonically 
increases with redshift, supporting the biased-galaxy for- 
mation scenario. 

3. We make sub-samples of LBGs at z = 4 binned ac- 
cording to magnitude, and calculate the correlation length 
for each sub-sample. We find that a brighter sub-sample 
has a larger correlation length (or larger galaxy-dark mat- 
ter bias): the brightest sub-sample (M ~ M* — 1) has 
r = 7. 91^7^00 Mpc while the faintest sub-sample (M ~ 
M* + 0.5) has r = 2A±Q 2 5 h^ Mpc. Similarly, we make 
subsamples of LBGs at z = 4 according to E(B — V) es- 
timated from the UV-continuum color, i' — z' . We find 



that the dustiest LBGs in our sample (E(B — V) ~ 0.4) 
have r = &.2~^ 4 '^hi 00 Mpc, which is marginally larger than 
the average. These dependences of ro on UV luminosity 
and dust extinction seen in LBGs are in accord with the 
fact that SCUBA sources, which are thought to be very 
dusty starburst galaxies, have a large correlation length, 
r ~ 10/ifoo Mpc. 

(b) Theoretical implications: 

4. We investigate properties of dark halos hosting LBGs 
based on the luminosity function (SDS V) and correla- 
tion function (section 3) of our LBGs, using the Cold 
Dark Matter model (Sheth & Tormen 1999; Sheth, Mo, 
& Tormen 2001). We estimate the mass of dark halos 
hosting LBGs at z = 4 from the the galaxy-dark mat- 
ter bias at a large separation, b g , and the number den- 
sity of galaxies from the luminosity function. We find that 
LBGs at z = 4 reside in dark halos whose mass ranges over 
10 11 - W 13 h^M Q , and that bright LBGs (M ~ M* — 1) 
have more massive dark halos (4.6l 3 '4 x ^ 12 h^o Mq) than 
faint LBGs (M ~ M* + 0.5; 1.7±o$ x IO'^^Mq). We 
compare the number density of LBGs with that of all dark 
halos existing at z = 4 (Figure 9) to find that the number 
density of massive dark halos (~ 10 13 h^M Q ) is compa- 
rable to that of LBGs while the number density of less 
massive dark halos (< 10 12 /i 7 " 1 M Q ) is 3 to 10 times higher 
than that of LBGs. This indicates that the majority of 
the low-mass (< 1O 12 /i 7o 1 M ) dark halos have no LBGs, 
and that most of the galaxies at z = 4 do not have UV- 
continuum emission originated from star-formation activ- 
ity which is strong enough to exceed our detection limit 
(M1700 ~ —20 or <~ 20/i 70 2 M Q yr~ 1 for extinction-corrected 
star- formation rate). Thus, we propose that LBGs are 
high-z galaxies which have episodic star formation and are 
detected only when star-formation activity is high. 

5. We perform a similar analysis to LBGs at z = 5 in 
our data and LAEs at z — 5 obtained by SDS II. We find 
that both LBGs and LAEs at z = 5 have dark halos with 
<~ 10 12 Ii^qMq. The number density of LBGs is about 2 
times smaller than that of dark halos with <~ 10 12 HjqMq. 
On the other hand, the number density of LAEs is about 
3 times higher than that of dark halos. Thus, dark halos 
of — 1O 12 /i 7 " o 1 M preferably host multiple LAEs. 

6. We apply a similar analysis to LBGs, X-band selected 
galaxies (FIRES galaxies), and SCUBA sources at z = 3 
whose correlation lengths and number densities are given 
by Giavalisco & Dickinson (2001), Daddi et al. (2003), and 
Webb et al. (2003), respectively. We find that LBGs with 
L > L* and FIRES galaxies have dark halos whose mass 
is ~ lO 12 ft- 7 " o 1 M , while dark halos for SCUBA sources are 
more massive (<~ 3 x 10 13 /i 7 " 1 M Q ), suggesting that SCUBA 
sources could reside in proto-cluster regions. The number 
density of LBGs is three times lower than that of dark ha- 
los existing at z = 3, while the number density of FIRES 
galaxies is about 4 times higher than that of dark halos. 
This implies that only about 1/10 (~ 0.3/4) of galaxies 
residing in dark halos with <~ 10 12 /i 70 1 M Q have an active- 
star formation whose star-formation rate (SFR) exceeds 
~ 5/ito 2 M Q yr _1 (~ 20/i 70 2 MQyr~ 1 if dust extinction is 
corrected), so that they are identified as LBGs by their 
bright UV continuum. This difference between LBGs and 
FIRES galaxies may support the idea that LBGs are galax- 
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ies which happen to be bright at UV wavelength due to 
episodic star formation. 

7. We find that the typical mass of dark halos hosting 
L > L* LBGs is about 1 x 10 12 h^ M©, which is com- 
parable to the total mass of the Milky Way Galaxy at 
the present epoch, and that the typical mass is almost 
constant over z = 3 to 5. It implies that the mass-to- 
luminosity ratio of LBGs (the mass of a hosting halo di- 
vided by the UV-luminosity of a LBG in the halo) does 
not change largely from z = 3 to 5 or, equivalently, the 
star formation efficiency is almost constant over z — 3 — 5. 

8. Using the CDM model, we estimate the mass of 
present-day dark halos hosting descendants of high-z 
(z = 3 — 5) galaxies. We find that the mass of the de- 
scendants ranges from lO 13 /if o 1 M to 10 15 h^M Q , which 
are comparable to the mass of present-day clusters and 
groups. Thus, dark halos hosting the observed high-z 
galaxies will evolve into clusters and groups. In other 



13 

words, most of the descendants of high-z galaxies studied 
here are probably member galaxies in clusters and groups 
today. The number densities of FIRES galaxies, LAEs, 
and faint LBGs exceed that of present-day galaxies with 
Mb g < —19. If most of these high-z galaxies are progeni- 
tors of present-day galaxies with M& < —19, they should 
experience mergers at least a few times from z = 3 — 5 to 
z = 0. 
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or 

http://hikari.astron.s.u-tokyo.ac.]p/ A ouchi/work/astroph/SDS_V_VI/SDS_Vl.ps 



Fig. 1. — Rcdshift distributions, N(z), of our LBG samples obtained in SDS V. The solid and dashed lines indicate the number-weighted 
completeness of the LBG samples in the SDF and the SXDF, respectively. Since the tighter selection criteria (eq. 4) are applied to Viz-LBGs 
of the SXDF, their selection window is narrower than that of the Viz-LBGs in the SDF. 
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or 
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Fig. 2. — Sky distribution of z ~ 4 Bffi-LBGs detected in the SDF. The filled circles indicate Bffi-LBGs; larger symbols mean brighter 
magnitudes (2"</> aperture magnitudes in i' band), as defined in the figure. The thick solid line outlines the area of B, V, i', and z' images 
after regions of low S/N are trimmed. The blank area in the lower-right corner outlined by the dotted line does not have ii-band data and 
thus the BRi-LBG selection is not applied. The shaded regions are large masked areas (> 2000 arcsec 2 ) where we do not attempt the object 
detection because of the presence of bright stars. The small panel inserted in the lower right corner presents a reduction of the main panel 
on the scale of one to five. The ten boxes drawn by dotted lines indicate the positions of ten CCDs of the Suprime-Cam. The real size of 
individual CCDs is 6. '7 X 13/3. The projected scale of 10/iJ~ 00 Mpc at z = 4.0 is shown in the bottom left of the panel. North is up and east 
is to the left in this image. 
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Fig. 3. — Sky distributions of z ~ 5 Viz- and ffiz-LBGs detected in the SDF. The black filled (red open) circles present Viz-LBGs 
(ffiz-LBGs), and the size of the circles increases with the apparent brightness in the z 1 band, as defined in the figure. The positions of Lyman 
a emitters (LAEs) at z = 4.86 ± 0.03 obtained by SDS II arc also plotted by blue crosses. The thick solid line is the whole field of view 
(FoV) of our observations in the SDF, and the dashed line indicates the border of the area of the ij-band image. The shaded regions arc large 
masked areas (> 2000 arcsec ) where we do not attempt the object detection because of the presence of bright stars. The projected scale of 
10/i^QQ Mpc at z = 4.7 is shown in the bottom left of the panel. North is up and cast is to the left in this image. 
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or 

http://hikari.astron.s.u-tokyo.ac.]p/ A ouchi/work/astroph/SDS_V_VI/SDS_Vl.ps 



Fig. 4. — Same as Figure 2, but for BRi-LBGs detected in the SXDF. North is up and east is to the left in this image. The small panel 
inserted is a reduction of the main panel on the scale of one to five. The positions of nine CCDs are shown by dotted lines. 
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Fig. 5. — Same as Figure 4, but for Viz- and Riz-LBGs detected in the SXDF. Note that a search for LAEs has not been carried out in 
the SXDF. North is up and east is to the left in this image. 
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Fig. 6. — Angular correlation functions, ui(B), of BRi-LRGs (top), Viz-LBGs (middle), and ffiz-LBGs (bottom). The filled circles are for 
the SDF, and the open circles in the top panel are for the SXDF. No significant signal is detected for either Viz-LBGs or Riz-LBGs in the 
SXDF. The solid lines show the best-fit power law function, u>{6) = J 4Jj iw #~' 3 , where is fixed at 0.8. Solid line of the top panel is the best-fit 
function both for the SDF and the SXDF data. 
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Fig. 7. — Correlation length tq (top panel) and galaxy-dark matter bias b g (bottom panel) as a function of redshift. The filled circles 
indicate our LBGs at z = 4 — 5, and the filled diamond is for LBGs at z = 3 (Giavalisco & Dickinson 2001; results based on the SPEC sample 
with R < 25.0). These LBGs at z = 3 — 5 have a similar luminosity (> L"). The open circle indicates LAEs at z = 4.9 obtained by SDS 
II. The filled pentagons, triangles, and squares arc tq of galaxies at z = — 1, derived by Loveday, Maddox, Efstathiou, & Peterson (1995), 
Carlbcrg et al. (2000), and Brunncr, Szalay, & Connolly (2000). The bias values of these galaxies at z = — 1 arc calculated in the same 
manner as given in section 3.4. The dotted line shows the correlation length of the underlying dark matter obtained with a non-linear theory 
(Peacock & Dodds 1996). The dashed line indicates tq of dark matter predicted by linear theory which we use in section 3.4. The solid lines 
indicate tq and b g of galaxies with B < — 19 + 5 log h predicted by the semi-analytic model of Kauffmann, Colbcrg, Diaferio, & White (1999). 
Our results for LBGs at z = 3 — 5 also show good agreement with the predictions of the semi-analytic model for galaxies with mjj < 25.5 
(Baugh ct al. 1999). 
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Fig. 8. — Magnitude dependence of the correlation length, ro, and the galaxy-dark matter bias, b g , for LBGs at z = 4 in the SDF. The 
four filled circles correspond to the four subsamplcs. The solid line shows the scaled Norberg's law. The ticks in the top margin indicate the 
number density of observed LBGs calculated from the luminosity function given by SDS V. The ticks in the right margin indicate the mass 
of the hosting dark halos calculated by the formula given by Sheth & Tormen (1999) using b g we obtain (see text). 
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Fig. 9. — Mass function of dark halos hosting LBGs at z = 4. The ordinate shows either the number density of LBGs in a given dark-halo 
mass range (filled circles), or the number density of dark halos (thin solid line). AM shown in the ordinate axis means the number density in 
the mass range of 1 magnitude, i.e., AM = M/ (2.5 log(e)) = 0.92M. The thin solid line indicates the mass function of all dark halos at z = 4 
calculated using the formula given by Sheth & Tormen (1999). The four filled circles correspond to the four data points shown in Figure 8. 
The thick solid line corresponds to a fit of the scaled Norbcrg's law to the four data points of Figure 8. 



24 Ouchi et al. 



See http:/ /hikori.ostron.s.u-tokyo.oc.]p/ A ouchi/work/ostroph/SDS_V_VI/SDS_Vl.pdf 



or 



http:/ /hikari.astron.s.u-tokyo.ac.]p/ A ouchi/work/astroph/SDS_V_VI/SDS_Vl.ps 



Fig. 10. — Same as Figure 9, but for galaxies at z = 3. The two filled circles show LBGs at z = 3 in the SPEC sample (right) and the 
PHOT sample (left) of Giavalisco & Dickinson (2001). The open circle with an arrow indicates the measurement for the HDF-N sample 
of Giavalisco & Dickinson (2001). The filled triangle corresponds to all galaxies found in FIRES (Daddi ct al. 2003) which is a ii"-band 
limited high-z galaxy survey with K < 24 (Labbe et al. 2003), while the two open triangles are for FIRES galaxies with blue color (left; 
J — K < 1.7) and red color (right: J — K < 1.7), respectively. The star corresponds to SCUBA sources whose correlation length is reported 
to be r = 12.8 ± 4.5/iT" 1 , Mpc by Webb et al. (2003) (see Scott ct al. 2002 for their number density). 
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Fig. 11. — Same as Figure 9, but for LBGs and LAEs at z = 5. The filled circle and filled square indicate LBGs and LAEs, respectively. 
The open square is also for LAEs, but when a higher dark-halo mass, 3.5 X lO 12 /^ 1 Mq, is adopted. See footnote of section 4.2.3 for details. 
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Fig. 12. — Same as Figure 9, but for present-day descendants of high-z galaxies at z = 3, 4, and 5. We calculate the dark halo mass of 
descendants using Sheth, Mo, & Tormen (2001), but assign to the descendants the same number densities as measured for galaxies at high 
rcdshifts. The blue, black, and pink filled circles show descendants of LBGs at z = 3, 4, 5, respectively. The red square is for descendants 
of LAEs at z = 5. The blue triangle and star mark denote all ii"-selected galaxies with K < 24 (FIRES) and submillimcter-selected galaxies 
with > 3mjy (SCUBA). The thick solid line corresponds to the scaled Norberg's law at z = 4 (Figure 9) evolved to z = assuming the model 
of Sheth, Mo, & Tormen (2001). The dotted line indicates the relation for 2dFGRS galaxies with M\ >g < —19, which is calculated from the 
mass function of dark halos (solid line) and the occupation number derived by van den Bosch, Yang, & Mo (2003) 
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Table 1 



Photometric Samples of Galaxies 



Field Name 


Sample Name 


Detection Band 


Number 


Magnitude Limit+ 


SDF 


Bffi-LBG 


i' 


1,438 


i' < 26.3 


SDF 


Viz-LBG 


z 


246 


z < 25.8 


SDF 


ffiz-LBG 


z 


68 


z < 25.8 


SXDF 


BRi-LBG 


i' 


732 


i' < 25.8 


SXDF 


Viz-LBG 


z' 


34+ 


i' < 25.3 


SXDF 


ffiz-LBG 


z' 


38 


i' < 25.3 



* Total magnitudes 

f V^-LBGs of SXDF arc selected by eq. (4). 
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Table 2 

Objects selected by two criteria 



Field Name 


Criteria^ 


Nobs" 


AW 1 


SDF 


Viz+Riz 


37 


41 ± 8.2 


SDF 


Viz+LAE* 


6 


4 ± 0.4 


SDF 


ittz+LAE* 


2 


2 ± 0.3 


SXDF 


Viz+Riz 


3 


5 ± 3.1 



'Two selection criteria. 

^Number of objects identified by both selection 
criteria. 

^Number of objects estimated from the contam- 
ination, completeness, and redshift distribution of 
the samples as described in the text. 

^"LAE means the criteria for selection of Lyman 
a emitters (LAEs) at z = 4.86 ± 0.03 in the SDF. 
Details of the LAE sample arc described in SDS II. 



Lyman Break Galaxies at z ~ 4 and 5 



Table 3 

Summary of the clustering measurements. 



Sample Name 


<2> 


Criteria* 




/<=(%)* 


A'™ 




BRi-LBG 


4.0 


all (»' < 26) 


1438+732 


2 0. 


qn +0.11 „ ,, + 0.9 
yu — 0.07 z -°-0.6 


2 4+ - 9 


BRi-LBG 


4.0 


all (»' < 26) 


1438+732 


2 


0.8(fix) 1.7+^i 




Viz-LBG 


4.7 


all (z 1 < 25. 8) 1 


246 


26 


0.8(fix) 2.1+^% 




Riz-LBG 


4.8 


all (a' < 25.8) 


68 


40 


0.8(fix) 0.4±| g 


1 r>+8.9 
1 - z -8.4 


BRi-LBG 


4.0 


i' < 25.3 1 


357 


1 


0.8(fix) 2.6+o g 


2 6+ UU 


BRi-LBG 


4.0 


23.3 < »' < 24.8 


103 


1 


0.8(fix) 6.O+3 2 


6 1 + 3 3 

— 3.2 


BRi-LBG 


4.0 


23.8 < *' < 25.3 


353 


1 


0.8(fix) 2.4±g;| 


2 4 +i.o 


BRi-LBG 


4.0 


24.3 <i'< 25.8 


724 


1 


0.8(fix) 1-3±°"| 


1 q+0.5 


BRi-LBG 


4.0 


24.8 < t' < 26.3 


1335 


1 


0.8(fix) l.2+_l% 


1 q+0.3 
1 -°-0.2 


BRi-LBG 


4.0 


-0.2 < E(B - V) < 0.2' 


253 


1 


0.8(fix) 2.9+J g 


2 9+ 1 3 


BRi-LBG 


4.0 


-0.1 < E(B - V) < 0.3' 


268 


1 


0.8(fix) 2. 0+12 


2-o+il 


BRi-LBG 


4.0 


0.0 < E(B - V) < 0.4 a 


228 


1 


0.8(fix) 3.2±^ | 


o n+1.5 
J - Z -1.5 


BRi-LBG 


4.0 


0.1 < E(B - V) < 0.5" 


153 


1 


0.8(fix) 2.6±1;^ 


o 7+2.2 
^■'-2.2 


BRi-LBG 


4.0 


0.2 < E(B - V) < 0.6" 


76 


1 


0.8(fix) 6.0±*'l 


6 1+ 4 ' 4 



t Selection criteria of the sub-sample, "all" means the whole sample of LBGs. 

Number of the sample galaxies. 
* Fraction of contamination in the sample obtained in SDS V. 
^ Correlation amplitude, A w , corrected for contamination. 

™ Clustering measurements for LBGs with L > L* used to calculate ro and 6 g plotted 
in Figure 7. Corresponding ro and b g values arc summarized in Table 5. 

a These E(B — \ / )-dcpcndcnt subsamples are constructed from 357 bright LBGs with 
i < 25.3. We calculate E(B-V) from i - z' color (section 3.6.2). The E(B-V) range of 
each subsamplc corresponds to —0.18 < i' — z' < 0.16 (subsamplc of —0.2 < E(B — V) < 
0.2), -0.10 <%'-%'< 0.24 (-0.1 < E(B - V) < 0.3), -0.01 < i - z < 0.33 (0.0 < 
E(B - V) < 0.4), 0.07 < i' - z' < 0.41 (0.1 < E(B -V) < 0.5). and 0.16 < i' - z' < 0.50 
(0.2 < E(B - V) < 0.6). 



Note. — Clustering amplitudes, A™ w and A w , arc in units of arcsec^. We regard A w 
as the best estimates. 
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Table 4 

Summary of the correlation length and the galaxy-dark matter bias. 



Sample Name 




Criteria^ 


r raw 










i> 9 


BRi-LBG 


4.0 


all (i' 


< 26) 


4 0+° 


2 

.2 




1 
1 


4.il° : 2 2 


2 9+° 1 


Viz-LBG 


4.7 


all (z' < 


' 25.8) 1 


4.2+; 


.0 
.2 


3.4i° 


7 

.9 




4 6+°' 9 


BRi-LBG 


4.0 


i' < 25. 3 1 


b.i+; 


(J 

.1 


3.5± U 


6 
7 


s.ii;- 


o r + U.t) 
oo -0.7 


BRi-LBG 


4.0 


23.3 < i' 


< 24.8 


7.9+ 2 


1 

.7 


5.3+; 


3 
7 


7 9+ 2 1 


5.3li;? 


BRi-LBG 


4.0 


23.8 < i' 


< 25.3 


4.7+J 



.2 


3 3+° 


6 
.7 


4.8+;° 


3.4+.°,;° 


BRi-LBG 


4.0 


24.3 < i' 


< 25.8 


3.4t°„ 


6 
.7 


2 5+°' 


4 


3 5+° :7 


2 5+°' 4 

^ - -0.5 


BRi-LBG 


4.0 


24.8 < i' 


< 26.3 


3.3t°„ 


4 
.4 


2 4+°' 


2 
2 


3 4+°' 4 


2 4+°' 2 


BRi-LBG 


4.0 


-0.2 < E(B 


- V) < 0.2 a 


5.4+i 


y 

.6 


3.7+ 1 ;- 


S 



5 5+ la 


3 8+"- i< 


BRi-LBG 


4.0 


-0.1 < E(B 


- V) < 0.3 a 


4.4+; 


.3 
.8 


3.i+; 


8 
2 


4.5t; : 4 8 


o 1+0.8 
dl -1.2 


BRi-LBG 


4.0 


0.0 < E(B - 


- V) < 0.4 a 


5.7+; 


3 

.7 


3.9+; 


8 



5.811- 


3.9+°/? 


BRi-LBG 


4.0 


0.1 < E(B - 


- V) < 0.5 a 


5-1 + 3 


.0 
.2 


3.5+ + 


2 
1 




3.6t + ? 


BRi-LBG 


4.0 


0.2 < E(B - 


- V) < 0.6 a 


8.lt 2 4 


8 

.2 


5.4+.+ 


7 

6 


8.2+.^ 





t Selection criteria of the sub-sample, "all" means the whole sample of LBGs. 

^ Clustering measurements for LBGs with L > L* plotted in Figure 7. These are repeatedly 
presented in Table 5. 

a These E(B — V)-dcpcndcnt subsamplcs are constructed from 357 bright LBGs with i' < 25.3. 
Wc calculate E(B — V) from i' — z' color (section 3.6.2). The E(B — V) range of each subsamplc 
corresponds to -0.18 <%'-%' < 0.16 (subsamplc of -0.2 < E(B - V) < 0.2), -0.10 < 
% - z < 0.24 (-0.1 < E(B - V) < 0.3), -0.01 < i - z < 0.33 (0.0 < E(B - V) < 0.4), 
0.07 < i' -z' < 0.41 (0.1 < E(B-V) < 0.5), and 0.16 < i' - z 1 < 0.50 (0.2 < E(B-V) < 0.6). 

Note. — ro and b g are the correlation length and bias parameter, respectively, calculated from 
the contamination-corrected clustering amplitude, A^, while rQ aw and b^ aw are those calculated 
from the raw clustering amplitude, A^ aw (see Table 3). We regard ro and b g as for the best 
estimates, ro and r£ aw arc in units of h^ QQ Mpc. b g and 6^ aw are defined at 8/i^ 00 Mpc. 
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Table 5 

Summary of the clustering measurements for LBGs with L > L* and all LAEs plotted in Figure 7. 



Sample Name 


<2> 






Criteria^ 








6/ 


(7„Gfl-LBG(SPEC)t 


3.0 




< 


-20.5 ~ M* 


(fl < 


25.0) 


5 0+ - 7 




BRi-LBG 


4.0 


-^1700 


< 


-20.8 ~ M* 


(*' < 


25.3) 


5 1 + 10 


o c + 0.6 
°-°-0.7 


Viz-LBG 


4.7 


Ml 700 


< - 


-20.5 ~ M* 


(z' < 


25.8) 


5 9+ 1 ' 3 




LAE" 


4.9 


Mi 


700 


< -19.0 (z' 


< 27. 


3) 


6 2+ ub 

D - Z -0.5 


4 9+ u4 



' ro and 6 g for BRi-LRGs, Vzz-LBGs, and LAEs arc contamination-corrected values. 

* ro and 7 of J7-band drop-out LBGs arc obtained from the SPEC sample of Giavalisco 
& Dickinson (2001). 

^ ro and b g arc obtained by SDS II. These are contamination-corrected values. 

^ These criteria correspond to L > L* for LBGs at each redshift (see SDS V). 



Note. — Values for £>i?i-LBGs and V iz-LRGs arc the same as those presented in Table 

4. 
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Table 6 
Mass of hosting dark halos. 

Sample Name (z) Selection Criteria (M)[h~gM e ] log(n) N oc<lup (M) (z = 0) [h~^ M ] AT min (merge) 



U n GR-LBG(L>BOT) t 


3.0 






R < 25.5 




1 4+ 26 


X 


10 11 


-3.4 ± 0.06 


0.04 


2.8±i;l 


X 


10 13 


1.4 


£/ n &R-LBG(SPEC)* 


3.0 






R < 25.0 1 




1 3+ 10 


X 


10 12 


-3.6 ± 0.06 


0.3 




X 


10 13 


1.1 


(7By-LBG(HDF) t 


2.6 






V 60 6 < 27.0 




< 1 X 


10 10 


-2.7 ± 0.06 




1 o+°- 3 








FIRES- All 


3.0 






if < 24 




9 9 +13 


X 


10 12 


-2.9 ± 0.03 


3.5 


X 


10 14 


6.6 


FIRES-Red 


3.0 


K 


< 


24 & J - K > 1.7 


8 9+ 3 ' 7 


X 


10 12 


-3.2 ± 0.05 


16.0 


2 +0 ' 6 


X 


10 14 


4.7 


FIRES-Blue 


3.0 


K 


< 


24 & J - if < 1.7 


3 5+ 14 


X 


10 11 


-3.1 ± 0.02 


0.3 




X 


10 13 


3.4 


SCUBA 


~ 3 






S s50 > 3mJy 




2-s±?:| 


X 


10 13 


-5.0 ± 0.5 


3.5 




X 


10 14 


0.1 


Bfli-LBG 


4.0 






*' £ 26 






X 


10" 


-3.2 ± 0.03 


0.3 


6.5+ u * 


X 


10 13 


2.8 


BRi-LBG 


4.0 






i' < 25. 3 11 




1 +0 ' 9 


X 


10 12 


-3.8 ± 0.05 


0.3 


9.6+fi 


X 


10 13 


0.9 


BRi-LBG 


4.0 


23.3 


< 


i' < 24.8 (i' = 


i 24.0) 4.6± 4 ;° 


X 


10 12 


-4.3 ± 0.1 


1.3 


2 1 + 12 

Z -!-l.l 


X 


10 14 


0.4 


BRi-LBG 


4.0 


23.8 


< 


i' < 25.3 (i' = 


; 24.5) &.6±\% 


X 


10 11 


-3.8 ± 0.09 


0.3 




X 


10 13 


0.9 


BRi-LBG 


4.0 


24.3 


< 


i' < 25.8 (i' = 


i 25.0) 2.ltj'| 


X 


10 11 


-3.3 ± 0.08 


0.1 


5.o±i;l 


X 


10 13 


2.1 


BRi-LBG 


4.0 


24.8 


< 


i' < 26.3 (i' = 


; 25.5) 


X 


10 11 


-2.9 ± 0.07 


0.2 


4 6+ ' 1 

*' u -0.9 


X 


10 13 


4.8 


Viz-LBG 
LAE 


4.7 
4.9 




z < 25. 8 1 
JVS711 < 25.8 




8 9 +f! J 

°- a -6.4 
1 1 + 05 


x 10 11 
x 10 12t 


-3.8 ± 0.3 
-3.2 ± 0.3 


0.6 
3.2 f 


, ,+U.B 
1 -°-0.6 
1 4+0-2 

L -^-a.2 


x 10" 
x 10 14 


0.9 
3.8 



* t/-band drop-out LBG samples given by Giavalisco & Dickinson (2001). SPEC, PHOT, and HDF mean the spectroscopic 
sample (R < 25), photometric sample (R < 25.5), and photometric sample of HDF-N (R < 27), respectively. 

* Maximum mass and occupation number of LAEs arc 3.5^q ^ x 10 12 hy^ Mq and 35 (see footnote in section 4.2.3). 
^ Samples of LBGs with L > L* at z — 3, 4, and 5. 

Note. — {M)[hjQ Mq] and (M)(z — 0)[hjQ Mq] are the mass of dark halos hosting high-2; galaxies and the typical mass of 
their descendants at z — 0, respectively, log(n) is the number density of high-2; galaxies in units of mag -1 hj Mpc~ 3 . N OCC up 
is the occupation number, and 7V mm (merge) is the minimum number of mergers required from high-2: to z — 0, assuming that 
the high-.s galaxies evolve into present-day galaxies with M& < — 19 mag. 
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Table 7 

Mass of hosting dark halos for LBGs with L> L* and all LAEs. 



Sample Name (z) 


Selection Criteria 


<M>[/ t - 1 M ] 


log(n) 


N occup 


(M)(z = 0)[h 




ln (merge) 


C7„Gfl-LBG(SPEC)*3.0 
BRi-LBG 4.0 
Viz-LBG 4.7 


R < 25.0 
i' < 25.3 
z' < 25.8 


1.3+J-" x 10 12 
1.0+™ x 10 12 
8.9+| : » x 10 11 


-3.6 ± 0.06 
-3.8 ± 0.05 
-3.8 ± 0.3 


0.3 
0.3 
0.6 


7.it a 2 -» x 

9.6+_H X 
1 3+ ' 5 X 


TO 13 
10 13 
10 14 


1.1 
0.9 
0.9 


LAE 4.9 


ATS711 < 25.8 


l.l+i! b 2 x 10 i2 


-3.2 ± 0.3 


3.2 


1 4+ 02 x 


10" 


3.8 



J7-band drop-out LBGs given by Giavalisco Sz Dickinson (2001). 



Note. — (M)[hj} Mq] and {M)(z — 0)[hjQ Mq] are the mass of dark halos hosting high-z galaxies and the typical 
mass of their descendants at z — 0, respectively, log(n) is the number density of high- 2: galaxies in units of mag -1 
hj Mpc -3 . N occup is the occupation number, and 7V mm (merge) is the minimum number of mergers required from 
high-2 to z — 0, assuming that the high-z galaxies evolve into present-day galaxies with M(, g < — 19 mag. 



